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wry Abstract 

In 1985, Edward Nelson, who formulated the theory of stochastic mechanics, made an inter- 
esting remark on Bell's theorem. Nelson analysed the latter in the light of classical fields that 
behave randomly. He found that if a stochastic hidden variable theory fulfils certain conditions, 
the inequality of Bell can be violated. Moreover, Nelson was able to prove that this may happen 
without any instantaneous communication between the two spatially separated measurement sta- 
tions. Since Nelson's article got almost overlooked by physicists, we try to review his comments 
on Bell's theorem. We argue that a modification of stochastic mechanics published recently by 
Fritsche and Haugk can be extended to a theory which fulfils the requirements of Nelson's analysis. 
The article proceeds to derive the quantum mechanical formalism of spinning particles and the 
Pauli equation from this version of stochastic mechanics. Then, we investigate Bohm's version of 
the EPR experiment. Additionally, other setups, like entanglement swapping or time and position 
correlations, are shortly explained from the viewpoint of our local hidden- variable model. Finally, 
we mention that this theory could perhaps be relativistically extended. 
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1 Introduction 

In 1935, Einstein, Podolsky and Rosen (EPR) wrote an article in which they denied that quantum 
theory would be a complete theory of nature [1]. Around 1951, Bohm gave a more testable outline 
of the so-called "EPR paradox" [2,3]. He described a thought experiment with one source that ejects 
particles having opposite spin to two spatially separate Stern-Gerlach magnets of variable orientations, 
see Fig. [TJ Then, in 1964, Bell published a theorem about this paradox in the form of an inequality. 
It made clear that hidden variable theories fulfilling certain conditions would contradict quantum 
mechanics. Bell called these conditions "locally causal" and explained them in detail in [4]. 

In his first publication on that topic, Bell wrote: "If hidden parameters would be added to quantum 
mechanics, there must be a mechanism, whereby the setting of one measuring device can influence 
another spatially separated device" and the signal involved has to "propagate instantaneously" [5]. 
However, Bell's first contribution underwent several modifications. Over the years, more and more 
instructive proofs of his inequality were constructed by him and others. A collection of all of Bell's 
fundamental articles can be found in [6]. Finally in 1969, Clauser, Holt, Shimony and Home (CHSH) 
brought the inequality of Bell into a form suitable for experimental investigation [7] . 

In 1985, Nelson tried to analyse Bell's theorem with full mathematical rigour. The result of his 
analysis was that Bell's definition of "locally causal" , which is the starting point to derive Bell's 
inequality, could be divided into two separate conditions. Both were necessary to derive the inequality, 
but only one of them has to hold if a hidden parameter theory does not include instantaneous signals. 
Nelson published his result twice [8,9]. (However, a small correction was added later [10], in order to 
make Nelson's theorem compatible with Mermin's presentation [11].) Furthermore, there are studies 
from other mathematicians, for example [12], which give further insightful analysis of Nelson's work. 
Unfortunately, Nelson's articles on Bell's inequality got almost overlooked by physicists. The reason 
for this might be that Nelson found his own interpretation of quantum mechanics [13, 14] to be at 
variance with the requirements of his theorem for a model without instantaneous signalling effects. 

This article is organised as follows: In section 2, we review Nelson's contribution towards a math- 
ematically rigorous understanding of Bell's inequality. The physical implications of that contribution 
are discussed in section 3. A recent modification of Nelson's stochastic mechanics [15] is the topic of 
section 4. It is shown that this modification fulfils Nelson's requirements for a local hidden variable 
model which is compatible with quantum mechanics. In section 5, we extend this version of stochas- 
tic mechanics in order to include spinning particles described by the Pauli equation. Entangled spin 
states are discussed in section 6, where it is shown in detail how their correlations emerge without any 
instantaneous signalling effects. In section 7, we review some arguments from the literature against 
Nelson's thoughts on Bell's theorem. We come to the conclusion that these arguments are unfounded. 
The article closes with section 8, where we mention the possible use of our theory as a reasonable 
foundation of quantum mechanics in curved space-times. 

The sections 2 and 3 will need some understanding of mathematical probability theory. For a 
general introduction to this theory, see [16-19]. 
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Figure 1: Drawing of an EPR experiment with Stern-Gerlach 
magnets 1, 2 and axes /2, V. The black dots illustrate the 
particles moving to the detectors. Their spin direction with 
respect to the z axis is indicated by the small arrows around 
them. 



Figure 2: Space-time diagram of 
an EPR experiment. The dashed 
arrows display particle trajectories 



2 Nelson's analysis of Bell's theorem 

2.1 The setup of the EPR experiment in theory 

An EPR experiment consists of two measurement devices 1 and 2 which are space-like separated, as 
well as a particle source in the intersection of their past cones. The source ejects j pairs of particles, 
where j € N and 1 < j < S 6 N, to the two detectors. The particles in each pair have opposite spin 
and opposite momentum. At 1 and 2, the direction of the spins is measured with a Stern-Gerlach 
magnet. These magnets can be rotated perpendicular to the axis x of the particle trajectory around 
an arbitrary direction. The direction of the magnets is called /2 for detector 1 and v for detector 2 (see 

Figs. OQ ©. 
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are called "spin trajectories" and denote the time dependent spin values of the j-th particles that are 
sent to detector 1 and 2. The outcomes in the state spaces of </>jit g and 4>j2t s are elements of sets 
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Cl'j = {t 7 |}- The sets Clji and Clj2 belong to sample spaces of two probability spaces 
whose definitions will be given later. 

As parameter of the spin trajectories, the vector z corresponds to the axis around which the 
described particle rotates. If it rotates in the left direction with respect to z, we have 4>ji(2)t, =1\ and 
a rotation in the right direction is defined by 4>j\{2)tg ={• The value t is the time during which the 
particle flies. Both particles start at time t = at the source and arrive at a time t = T, where T > 0, 
in a Stern-Gerlach magnet. 

When the particles arrive in the Stern-Gerlach magnets at time t = T, the magnetic field of these 
devices could change the particle's properties, including the spin. In an EPR experiment, the parts of 
the detectors that might influence the particles are the axes \1 and v. Those axes can be chosen by the 
experimenter freely at will. Due to the Stern-Gerlach magnets, we have to describe our measurement 
results by device dependent random variables. For the two detectors 1 and 2, we define two axis 
dependent families of random variables Dj\p, '■ &ji >— ► fi'- 1-D and Dj2V '■ ^ 
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of these random variables are Q'- 
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Q,j2- The outcomes 
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Q!. = {t,4}, and for simplicity, we will assume that 
QjiD = Qji and flj2D = Qj2- Since the random variables Djip and Dj2V depend on the setting of the 
Stern-Gerlach magnets, we write their outcomes on the state space using the axis dependent notation 
ajifl and a^v- However, the way, how the events generated by Djip and D^v are connected to the 
events generated by 4>jit 3 and 4>j2T s will be postponed at the moment. 

The results at the detectors could depend on the preparation of the particles by the source. We 
say that the corresponding events happen at preparation stage. They are generated by the family of 
product random variables 



Sj\2 = <f>jl,t=0z <8 <f>j2,t=Q s ■ (fijl X %2) ^ (f^-i X ^2) • 



(1) 



The axis dependent family of product random variables 

XjijX2P = D jip (g) Sji2 <X> D j2 : (fljlD X Qji X VL j2 X Qj2D) !->• (fij-izi x fi^ x Oj. 2 x ^ 2Z) ) (2) 

generates events on an enlarged probability space. We define the sigma algebra T on that space by 
the power set T = V (£Ijid x Clji x Qj2 X £Ij2d)- The probability measure of the enlarged probability 
space will be denoted by P. 

Now, we introduce a new family of random variables 

4>jljt <X> 4>j2i? ■ (%1D X flji X Q j2 X Oj2r>) !-> (f^lD X ^2Z)) ) ( 3 ) 

whose definition is 

{-Xji/Ki? € (o-ji/r x O^ x Q,' j2 x a,^)} = {^1^ ® ^2? G (ctji/j x ct^,-)} ■ (4) 

It provides information about outcomes which happen simultaneously at the two detectors. The 
probabilities of the axis dependent events in Eq. (|1]) will be abbreviated by 

P ({<?W ® <t>j20 € (^-i/J X CTj 2i ;)}) = P^ijKS^ap (<^1/T X C^) . (5) 

The expression Pf r gf K defines an axis dependent family of probability measures for events of the 
form (Jjip x Uj2v- We say that these events are generated by the family of random variables fijip ® 4'j2v- 

There may be situations where we are only interested in the outcomes at a single detector. For 
this reason, we define the notations (Jjip, = <Jjip X £1'j 2 d an d a j^ = tym x ^j2V- The events (Jjip, =t 
= (crji/J =T) x Q'j2D an d Oji/J =1 = {vj\p, =4) x f2' 2 £> generate a sigma algebra which we denote by 
Fjip,. Similarly, the events a. } 2V =T = ®>jid x (&j2i> =T) an d cr^i? =| = n^ 1X) x (c^ =|) generate the 
sigma algebra J-j2u- From these definitions, it follows that Tj\p (or Tj2v) contains the information 
about the outcomes at detector 1 (or 2) only. 

Since the family of random variables 4'j\p®4'j2v is defined on the enlarged probability space, it will 
give us the possibility to investigate the relationship between the outcomes at the detectors and the ones 
at preparation stage. We let a third sigma algebra J-js contain the information about the outcomes at 
the source. Accordingly, T^s is generated by the events S\ = \Xj\p2v £ (fi' 1D X f x J. xri' 2D )} and 

s 2 = {x jlfi 2s e (a jlD x 4 x t xpf j2D )}. 

We can condition the events in T with respect to the sigma algebra J-js- The sigma algebra Tjs is 
generated by two disjoint events S\ and 5*2. The probability of each of these events to happen is > 
and their union is the sure event of the enlarged probability space. In this special case, the conditional 
probability of an event A £ J 7 with respect to Tjs is given by (see [16] at section 33, or [19] at page 
126): 

PM^^g^j^lxdPl*. (6) 

Eq. ((6]) is the definition of a random variable for outcomes in the enlarged probability space. It can be 
interpreted as the revised probability of an event A to happen with respect to the extra information 
about which events in Tjs occur for a given outcome. Obviously, P(A\Fjs) fulfils the following 
properties: 

• For all events iGf, the random variable P (A\ J-js) is measurable with respect to Fjs- 

• We have < P (A\ Tjs) < 1 for all events A G T and outcomes on the enlarged probability 
space. If A is the sure event Qjid x ilji x $7j2 x ftj2D, the conditional probability P (A\ Tjs) 
is unity for all outcomes on the enlarged probability space. In case A is the impossible event 0, 
the random variable P (A\ Fjs) is zero for all outcomes. 



• For every sequence of pairwise disjoint events (A n ) n&fi G T, the equality P ((J n =i -^"1 3~i8) = 

En=l P ( A n\FjS) holds. 

• Since the union of Si and S2 is the sure event, the unconditional probability of any event A ^ T 
can be computed with the expectation value EX [P ( A\ Tjs)] — J P (^4| Fjs) dP = P (^4). 

• For any event A £ T and another event Xjg £ J~jSi where Xjg 7^ 0, it follows from Eq. ([5]) that 
P ( A P| Xj?\ Tjs) = P (^4 I Fjs) for outcomes in Xjg, and P (j4p|Xj,?| Tjs) = for outcomes not 
in Xjg. 

In case that the events on the enlarged probability space are generated by <f>jip <8> 4'j2i7, we will use 
the following notation: 

P ^iA®*iaa (^lM X *j2ff\ -^s) = P ({ <Aj1/T <8> <^2i7 € {Zjljt X cr^tf)}! -^js) ■ (7) 

Now, all the mathematical structures needed to analyse the EPR experiment have been defined 
and we can proceed with the necessary locality conditions. 

2.2 Active locality 

At first, Nelson defined two different forms of locality: Active locality and passive locality. The meaning 
of active locality is: 

Whatever axes the experimenter selects at one measurement device, e. g. at 2, it does not change 
the probabilities of an event happening at 1, as long as 1 does not lie in the future cone of 2. Hence, 
active locality corresponds to a "no-signalling" condition. If it is violated, one can send instantaneous 
signals between two spatially separated locations. We can mathematically define active locality as 
follows: 

We let the area of measurement station 1 be disjoint from the future cone of station 2. With 
the axis /2 of 1 being left constant, different axes v 7^ V of Stern-Gerlach magnet 2 are chosen. The 
two pairs of axes influence the measurement results at the two devices through the random variables 
<Pjip <8> 4>j20 and 4>jiji ® 4>j2v' ■ We define a theory to be actively local if 4>jiii ® <t>j2P an d 4>jip, ® 4 > j2V 
are equivalent for probability of an event (Jjip, = Vjiji x Q'j2Di or 

P 4>jip<»4>j2a ( a jip) = P<t> j i(i<®4> j 2,7> ( a jlp) ■ ( 8 ) 

A similar condition should hold for the events Uj2v = ^' 7 id x &j2v at Stern-Gerlach magnet 2, provided 
that the spatial region of 2 is disjoint from the future cone of magnet 1: 

FVii*®^a» (°j'2«?) = P 3lfJ /«>0 3 2P (< 7 j2p) • (9) 

In Eq. ([9]), we have p! 7^ jj, for the axes of detector 1, with the axis v of magnet 2 being left constant. 
By Eq. (|7|), we can condition the probabilities of <jjip £ Fjifl and a 3 2v G 3~j2V with respect to 
Tjs- The events in TjS contain no information about the results at the detectors. Instead, they only 
define, what has happened at the source. The latter is located disjoint from the future cones of the 
two measurement devices. Therefore, in an actively local theory, the probabilities of the events in 
Tjs cannot be changed by the detector's axes. Accordingly, it follows from Eqs. ([6]) and (JT]) that 
p 3 ijr®</>j2s (vjiftl^js) and P^^fgi^aa Wj2p\ Fjs) wm have similar properties as in the equations IJHJ) 
and ©: 

' p <t>jiji®<t>w>{vji-il\Fis) = p <t> j i f t<8<t> 3 2,7' ( a jip\Fjs), (10) 

'P<t> 1 xp®<t>j2i>{ (X 32p\f3s) = P <t> :jli t>®<i> i 3e( a 'j2p\Fjs)- (11) 



2.3 Passive locality 

In this subsection, we consider the joint probability of (Jjipt G 3~jia & n d <Jj2v G ^2i7- It gives information 
about events which happen simultaneously at the spatially separated locations 1 and 2. We say that 
passive locality holds if 

P <W»<^2,7 [ a 3lP f] °3' 1V ■ F ^ S ) = P ^l f J«"/» J 2? (Vjlitl^js) P(t> j iP®<t>j2,7 ( a j2p\ FjS) , (12) 

for every pair of axes jl and v. 

It is possible to have active locality without passive locality. Furthermore, a theory that violates 
passive locality does not have to incorporate any signalling mechanism between the spatially separated 
measurement stations 1 and 2. However, the true physical implications of the passive locality condition 
will be postponed at the moment. 

2.4 Nelson's theorem 

With the definitions above, quantum mechanics fulfils the following properties: 

P 4w»te [Vjlp =t Pl ^ = -U + P £fi/r®£f2/z ( a jlfl =1 f] a 32p =t) = 1, (13) 



|E (p, V) + E (p, V) + E [pi, ff)-E (/i', &)\>2. (14) 



Equation (fl3j> implies that if the axes of the Stern-Gerlach magnets are the same, the spin values 
measured at 1 and 2 are always opposite. Using P^ 1<r ®^ 23 {crjip =T) + "^ 4>jip®4> 3 2p (°ji/J =1) = 1 an d 
P^ lfJ <»0 32f r ((Tj2p, =T) + p 3 i f r«i032 f j (°j2/z =1) = 1, we can rewrite Eq. ([13]) in the following form: 

p 3 ifj«><te ( CT ii/i =T) = P4> j ip®<P J 2p ( a jip =T Pi ^'2/2 =4) = P<t>jip®4>j2p ( a j2p =i) > 

+ + + (15) 

P<t>jip®<l>j2p ( a jlP =i) = P<l>3ip®4>32P ( a j^P =i H ^'2/2 =T) = P<t>jip®<t>j2P ( a j2P = T) ■ 

This makes clear that the events at the detectors are equivalent if the same axes at the separated 
measurement devices were chosen. 

The function E(p, v) in Eq. (fl4|) is called correlation coefficient. It is defined through: 

E (p, v) = P^®^ (ajip =T P| <?j20 =TJ + p J i A «.0 j2 p (°ji0 =4 n a ^ i) = y 

- P<t>jip®<t>j20 [Vjip =T H ^ 2 ' 7 = J-J ~ P 4>mt9<t>j2P \ a nP =1 f] a i2P =TJ ■ (16) 

It was shown by Nelson that Eq. (|14[) is in conflict with theories where Eq. (fT"3"|) and active as well as 
passive locality hold. Nelson's proof [8,9] (with corrections in [10]) proceeds as follows: 

We are interested in the properties of the conditional probabilities of the events <7jip € 3~jip, and 
<Jj2v S Fj2v given the sigma algebra J-js- At first, we will investigate, what the equivalence property 
of Eq. (|15[) combined with the assumption of passive locality implies for these conditional probabilities. 

If the axes at the two detectors are the same, i. e. p = P, we have, due to passive locality: 

p 0iip«><te \Vjip\\vj2p, Fjs) = P<t> n p®<pj2p ( a jip\^ F js)P<t> jl p<»4> j 2p Wj2p,\J : js)- (17) 

Since < P^®^ {cr j2 p\ Fjs) < 1, we get 



"<t>jlp®<t>j2p 



\ a jip n cr J 2 ' :r ^j - ^hi^^p (<rjip\fjs) ■ (is) 



With the definitions of <Xji# and (Jj2p, it follows from Eq. ([5]), the properties of the conditional proba- 
bilities stated in section 2.1, and Eq. ([7]) that 

P ^i*®^23 foi/*) = P ({<?W ® </W € CTji^}) 

= EX [P ( {^ <g> i2jS G crji^}! ^,-s)] 
= EX [P^ I? ®^ 33 (o-ji/jl ^-s)] 
and, similarly, 



^ 



:;.? 



(19) 



(20) 



Using Eq. (|15[1 . we observe that an event (Tjip, at detector 1 implies an equivalent event <Jj2p. at 2, i. e. 

Hence, the expectation values in Eqs. I|19p and (I20[) are equal. Accordingly, we get with Eq. (JT5J): 

p <t> 3 ip®<l>j2ii y<7jiil\\<7j2jt FjS) =P<t> J xp®<l> j2 p((7jiii\3 r js), (22) 



and 



3 (^ImPI ^ "^J - p ?W®te (c^l-^s)- 



(23) 



Now, we will look at the consequences of active and passive locality for the conditional probabilities. 
Passive locality demands for events <Jjia =t at 1 and cr^ =1 at 2 that: 



P J lM«"/'.,2P ( Vjlp =t f]^' 2 ' 7 _ ^ 



•^sj = P<t> jip <s»i> j 2i,( a jiii =Tl^s) 

xP^u,«^»(o-ia?=i|^s). (24) 

Active locality implies that the actions happening at 2 can only affect anything in the future cone 
of 2. Since 1 does not lie in that cone, nothing what happens in 2 can affect the events measured at 
1. Therefore, the random variables 4>jip ® 4>j2u and 4>jip. ® 4>j2p. are equivalent for the probability of 

P ^-i?®^2s ( a jip, =t\^js) = P<t> jiri ®<p j2 p (vjifl =T| Fjs) = p fl- (25) 

An analogous expression is true for the events at 2: 

P«W®0j2* {<7j2i> = 11 FjS) = P^i*®^ 2£ , (cTj2,7 =4| .Fjs) ■ (26) 

Without loss of generality, we can select an axis jl = v with v from Eq. (]26f> in the Eqs. (|X5f>. ([22)) 
and (|23f . It then follows from Eq. (|15|) that the event at detector 2 in the right hand side of Eq. (|26| 
implies an equivalent event of <Jj\v =T at 1. Using the Eqs. I|22|) and Q23[) with an axis /2 = v, we can 
conclude that 



P <Pjifl®<Pj2,y ( a o2u =11 ^5) = P^ii7®^ 2 £, (<7j2i? =11 ^Ss) 

= PtjipQtop WjW =T| -^s) = -Pi7- 

Plugging Eq. ([25]) and Eq. ([27]) back to Eq. pi)), we arrive at 



P ' <j> np <gi<t> iM \VjXp. — T \) (7 j2P —I 



J" 



js 



1 fl- 1 v 



(27) 



(28) 



The events ctj2i7 =T and crj2i? =1 are disjoint and their union is the sure event. According to section 
2.1, the sum of the conditional probabilities with respect to a sigma algebra is equal to unity for such 
events. Therefore, we get with Eq. (T2~7|) : 



2 A°J2V=1\fjs) = l-Pu. 



(29) 



Since active locality holds, a similar result can be written for the event Cji M =| at 1: 

= l-P^®^(o-ii3=T|^is) = l-^. (30) 

Due to passive locality and the Eqs. ([30]) and ([29]) . we have for the events (Tjip, —[ and (jj2ff =T : 

P^®^» (*ii/r =1 f>^ = T | ^s) = (1 - P/:) (1 - P*) ■ (31) 

In the same way, we can derive the relations 

P^,-fe (^i/z =T fl °W =t| ^s) = Pfl (1 - Pi;) (32) 

and 

PtunQfae [vjifl =4 fl G ^ =4 ^ s ) = ( X ~ P t) p z- ( 33 ) 

Using the Eqs. ((221), (O, @ and l[3"I]). we may define the function 

E(^z7|^ s ) ee P /I (l_P !7 ) + (l_P / j)P !7 _P / jP ? _(l_P /I )(l_P i7 ). (34) 

The conditional probabilities in Eq. (CM)) are all in the range < P^ ri a®^ i2i7 (cryi/j H <7j2i/| -^js) < 1- 
Therefore, the following inequality can be computed with four arbitrary axes fl, jl' and £>, i/ at the two 
stations 1 and 2: 

|E (ft, v\ T ]S ) + Eift^l^js) + E(i?,0\ T jS ) - E {fl 1 1 f\T ]S )\ < 2. (35) 

If we replace the conditional probabilities in Eq. Q35p by their corresponding unconditional probability 
measures, we arrive at the correlation coefficient from Eq. (|14p . The unconditional probabilities are 
also in the range of < P^ ;i 2<g>0 !2P {cjt/j, f] <Jj2v) < 1, and they are given by the expectation values of 
the conditional probabilities. Hence, an inequality analogous to Eq. (|35[) must be true for them: 



|E {fl, u)+E {fl, V) + E {flf, u)-E {fl!,V)\ < 2. (36) 

Equation (f3"r?)) is called Clauser-Holt-Shimony-Horne inequality [7]. It is a version of the statement in 
Bell's first article [5] 

|E(/J,i7)-E(/?,i7)|-E(i?,*')<l, (37) 

which can be similarly derived. Both Eq. (|36[) and Eq. (|37[) are violated in quantum mechanics and 
this violation was confirmed experimentally in 1982 [20]. 

The analysis of Nelson leads to two possibilities: Since Bell's inequality can only be derived assuming 
that active and passive locality holds, in a theory that reproduces quantum mechanics, either active 
locality or passive locality has to be violated. A violation of active locality would mean that action 
at a distance is possible. The physical implications of a violation of passive locality will be explored 
in the next section. 

It should be remarked that Jarret arrived at a similar conclusion [21], just one year before Nelson 
published his articles. Although Jarret was not able to identify the correct probability space given 
by quantum mechanics for the EPR experiment, one can find a preliminary definition of active and 
passive locality in his publication. 



3 What the failure of passive locality implies — a deeper anal- 
ysis 

Stochastic processes can be constructed such that they act on each particle actively local, but still are 
violating passive locality. This was shown by Nelson with two examples [8,9]. The first example was 
a rather ad-hoc construction of a stochastic process fulfilling the Klein-Gordon equation. The second 
example, mentioned briefly at p. 537 in [9] and, with some more detail, at p. 450 in [8], was a field 
theoretic variant of stochastic mechanics. It was originally developed by Guerra and Ruggiero [22], and 
it allowed to interpret various bosonic field theories, including the Maxwell field, the vector meson field 
and the linearised gravitational field, as the result of a stochastic process [23-25] . Unfortunately, this 
theory had no tempting particle interpretation as Nelson's original theory had, where the particles move 
on stochastic paths. Therefore, the program of Guerra and Ruggiero was never developed extensively. 

In 2003, Fritsche and Haugk published another modification of stochastic mechanics [15]. Unlike 
the model of Guerra and Ruggiero, the work in [15] constructs a non-relativistic hidden variable theory 
for fermions. The present article aims to show that the theory of Fritsche and Haugk also violates 
Bell's inequality without any communication between the two separated measurement stations of the 
EPR experiment. However, in order to reach our goal, we need to evaluate the implications of Nelson's 
theorem further. 

We begin by analysing some deeper consequences of passive locality. Faris showed in [12] that if 
passive locality is combined with relation (|13[) from quantum mechanics, it immediately leads to another 
condition, which he calls "deterministic passive locality" . Suppose two spatially separated events 
happen or not happen always together with probability 1. Then, deterministic passive locality states 
that these events are determined by some random event at a prior preparation stage. Deterministic 
passive locality is defined mathematically as follows: 

Let the events (Jjifi <E Fjip, a t 1 an d crj2fi £ Tjiji at 2 be equivalent with respect to P^ .^^ . 2 , and 
passive locality hold. Then, there must be an event Xjg £ Tjs at preparation stage, which is equivalent 
to both ctjip, and <Jj2pi- 

In his contribution, Faris states that a similar result is presented by Redhead in [26] at pp. 101- 
102. Redhead claims, it would have been discovered at first by Suppes and Zanotti [27]. Recently, 
another derivation was given by Conway and Kochen [28] (who called their result "free will theorem" to 
emphasise its implication of the particle's random behaviour). They relate the statement given below 
to the famous Kochen-Specker theorem [29]. The articles of Conway, Kochen and Specker are using 
geometric arguments. In contrast, the proof of Faris relies on probability theory, and explicitly uses 
the property of passive locality. The two different methods to prove the same lemma demonstrate an 
interesting connection between Bell's theorem and the one of Kochen-Specker. That the inequalities 
of Bell and the theorem of Kochen-Specker are in some way related to each other was also pointed out 
by Mermin [30] on the grounds of a different argument. The derivation below will follow closely the 
lines of Faris: 

As shown in section 2.4, it results from the equivalence property of Eq. (fl~5|) and the passive locality 
condition of Eq. (TT2|) that 



^<t>jiji®<Pjip 



y (7 jip\\< 7 j2fl Fjs) =P<t> j ia®<P J 2p( a jifi\J 7 js) (38) 

and, similarly, 

p 0ji^«>fe \ a jip. [) a &it Pis) = P <W»«W {<Tj2jt\Fjs) • (39) 

Using the assumption of passive locality again, we get with Eq. ([38]) and Eq. (f39|) : 

= (Pfe^te(^ipl^s)) 2 . (40) 



Eq. (|40|) implies, that the random variable P^^-igxp^a {cjifil J~js) can only have the values 1 and 0. 
The union event of the outcomes in the enlarged probability space for which P^- ljz ®^- 2 - (cjiu\ J~js) bas 
the value 1 will be denoted by Xjg . From this definition and Eq. (|40p . it follows that P^ j-igi^ 2t (<7jiu\ Fjs) 
is equal to the indicator function of Xjg\ 

p 4w»<te ( a jip\ Fjs) = Ixjg- (41) 

According to section 2.1, P^,^.^ {&jifi\ Fjs) is measurable with respect to Tjs- Therefore, we must 
have Xjg e fjs- 

Now, we recall the definition 

p «W®te (vjiillfjs) = p ({<?W <8 fe/r G o-ji^l-^s) (42) 

from section 2.1. In Eq. (|19p . the unconditional probability of Oji£ is given by the expectation value 
p 0jiM«"te ( a jiii) = EX [P^iM«"^i2M ( a jip\ Fjs)] 

= ex [P ( {^ ® i2/ j e a il? }| ^-g)] 

= / P ( {<W <8> ^2/z G ^i/rll ^s) dP. (43) 

Similarly, the probability of an event Xj? on the enlarged probability space with its measure P can be 
computed by the expectation value of the indicator function l x .-,: 

EX [1 X .J = j l x „ dP = P {xjg) . (44) 

Using the Eqs. (43]), (42]), (41]) and g4]), we can conclude that 

P<te®<te {vjiit) = p (#j*) ■ (45) 

In section 2.1, we also have learnt that in case of Xjg ^ 0, the conditional probability of the intersection 
between the two events {4>jip <8> 4>j2p G CTji/j} and jcj^ € •Fjs is equal to 



p ({<?W ® <^2/2 G crjip} Pi a^-j /}. 



P ({<?W ® <Aj2# G CTji/i}! J^s) for outcomes in a;^, 



for outcomes not in Xjg, 



(46) 



On the other hand, we get with Xjg = 0: 

P f {4>jip <8> ^j2/j G o-jip:} P| xj-A Tjsj = P (0| .Fjs) = for all outcomes. (47) 

The Eqs. US]), (42]) and flU) lead to the relation 

p ( {</W ® <?W G crji^} P x J2 - Tjs) = P^ lp ^ m ( o-ji^l F jS ) = lx^- (48) 



Due to Eq. (47|) and the definition of Xjg, Eq. (|48|) also holds in case of Xjg = 0. Computing the 
expectation value from both sides of Eq. (|48|) yields 



p ({<?W ® <?W e (^1m x fi j-2c) } P &i*) = p (Xjz) ■ (49) 



Eqs. (|45[) and (|49[) imply that the event 0ji£ G -Fjiit is equivalent to an event Xjg <E •T'j's- It follows 
from Eq. (|38|). Eq. (39]) and the second line of Eq. gOJ) that 



^4>nit®^m ( CT iipl -^is) - (P^®^ (o-j2/r| fjs)) ■ (50) 
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Using Eq. (|4Tj) together with Eq. ([50]) . we arrive at 

p <W»te {&jip\ Fjs) = P<f>3n;®<i>32ii (°j'2#| Fjs) = lx i? - (51) 

A similar calculation as in the Eqs. (|43 |) - (j49|) but for Oj%p. then shows that dj2p € ^2/1 is equivalent 
to the same event Xjg € Tjs which ajip e -^i/j is equivalent to. 

Deterministic passive locality states that the events at the detectors must be predetermined by 
another event at preparation stage. The failure of this condition implies for a hidden variable descrip- 
tion of the EPR experiment that during the time of its flight from the source to the measurement 
apparatus, each entangled particle has to undergo a stochastic process which constantly changes its 
spin orientation from up to down and vice-versa. This must happen such that the particle's spin value 
at the detector is not determined by the spin value the particle has had at the source. 

Due to Eq. (|13p , the two spin values in each entangled particle pair must be opposite at every time 
until the particles arrive at the detectors. It may be that the experimenters had chosen the same axes 
of the Stern-Gerlach magnets at the two separated measurement stations. In this case, the spin values 
of the entangled particles must remain to be opposite after the particles have passed these magnets. 
Finally, for a theory compatible with relativity theory, all procedures involved must preserve active 
locality. 

If these conditions can be fulfilled, Bell's inequality will be violated without instantaneous signalling 
between the two measurement stations 1 and 2. As shall be shown in the next sections of this article, 
there is no reason to believe that hidden variable theories are unable to do so. 

4 Locality in the hidden variable theory of Fritsche and Haugk 

4.1 Derivation of the single-particle Schrodinger equation from Brownian 
motion 

In 2003, Fritsche and Haugk published a modification of Nelson's stochastic mechanics [15]. In contrast 
to Nelson's original contribution [13, 14], the model of Fritsche and Haugk is an actively local theory 
with a non-markovian stochastic process. Moreover, it can, with a slight modification, be made 
to violate passive locality. Here, we will not discuss in detail how the (many-particle) Schrodinger 
equation (and the corresponding quantum mechanical operators) can be derived from the stochastic 
process involved because this is fully described in the original article of Fritsche and Haugk. Instead, 
the focus of this section is on superposed states, and we show the theory to be actively local even 
if many-particle systems are considered. In turn, we will have enough material to explain how the 
correlations in entangled states emerge. 

The starting point to derive the modified stochastic mechanics of [15] is to interpret the energy 
time uncertainty relation AEAt sw H/2 as an indication that in quantum mechanics energy is only 
conserved strictly in physical interactions, but not otherwise. Particles are assumed to be under 
statistical energy fluctuations. In the theory of Fritsche and Haugk, there occurs a deviation AE of 
the particle's energy E during a time interval At, after which the initial energy E has to be restored 
on the average. In other words: Kinetic energy is thought to be conserved on the average only. The 
second assumption is in spirit of Ballentine's famous statistical ensemble interpretation [31]. In this 
interpretation, the quantum state vector \tp) does not represent an individual particle, but a statistical 
ensemble of infinitely many of them. 

For this reason, the probability to find the electron of a hydrogen atom in a volume £1 at time t, 
given by 

p(Q,t)= /V(r,ihMr,i)d 3 r, (52) 

Jn 

where ij)(r, t) is the electron's wave-function, describes the following situation: In order to get all the 
information of p(Cl,t) experimentally, one has to prepare infinitely many identical hydrogen atoms in 
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the same state, including the same orientation in space. Then, one has to measure at each atom if 
one finds the electron in £1 at time t after the preparation of the state. Finally, one must compute the 
electron's probability of being there with all the results from the infinitely large ensemble. 

Accordingly, one has to deal with an infinitely large number of sample trajectories in order to 
describe the probabilistic time evolution of the statistical ensemble of electrons described by a quantum 
mechanical system. The idea of energy conservation on the average for the individual particles of this 
ensemble is implemented in [15] as follows: Fritsche and Haugk divide all the sample trajectories into 
two sub-ensembles with equal average velocities, and therefore the same average kinetic energy. In 
the sub-ensemble {A}, each particle undergoes a Brownian motion with a negative friction coefficient 
during a time-interval At. The particles in the sub-ensemble {B} underlie the same process, but with 
a positive friction coefficient. After each time step, some of the particles are interchanged between 
{A} and {B}, until the kinetic energy of these sub-ensembles agrees again. At this point, the entire 
procedure is repeated. 

For a single particle of rest mass mo, we start our discussion with two Langevin equations for j-th 
trajectories: 



molfit) + ^xf(i) = Ff \t) + Ff °™(t), 
m x A (t) - —kf(t) = Ff *(£) + Ff own (t). 



(53) 



The three dimensional vector x- ' (t) = IxJ- '(4),x- 2 ' (t),x 3 (t)\ denotes the hidden position 

variable for the j-th sample trajectory of a particle in the {A} ({B}) ensemble at time t, where j E N 
and 1 < j < S € N. The time dependent trajectories are defined on a coarse grained timescale r 
at which the Brownian motion process takes place. The vector F° xt (i) = (F^ t (i),F°| t (i),F°3 t (t)) is 

an arbitrary external force field and F Brown (i) = (Ff{ own (t),Ff™ wn (t), F B 3 rown (i)) denotes the random 
force the j-th particle is subjected to. This random force is assumed to have a Gaussian distribution 
in each coordinate k. The probability of an outcome F^f own (£) is given by 



P (Ffr n (t)) = -^^f=e _ V v^b F -° Wn(t) J . (54) 

V 71 " V m okBT 



In Eq. (|54p . t co u is a mean time of momentum transfer, /cb is Boltzmann's constant and T is the 
temperature of the heat-bath that surrounds the j-th particle. The Gaussian distribution defined by 
Eq. (|5^|) should be independent for every coordinate 1 < k < 3 and trajectory j of the random force 
field. 

One can separate the sample trajectories x ■ into a sum caused by the random and the external 
force: 

xf B ) W= x^(t) + x^ (i ). (55) 

Since F^ rown (£) is independent in each coordinate, there does not exist any dependence between dif- 
ferent components of the velocities xi (t) , i r - 2 (t) and x r ^ 3 (t) , as well as there is no dependence 
between x.' \t) and x - l '(t). With this property in mind, Fritsche and Haugk average in {A} and 



? A(B) 



{B} about all j-th sample trajectories x ■ (t) after the two sub-ensembles were made infinitely large 
(S — > oo). Finally, they arrive at two differential equations: 

| (v B - u B ) + (v B + u B ) • V (v B - u B ) - „A (v B - u B ) = -Li?- 

| (v A + u A ) + (v A - u A ) • V (v A + u A ) + uA (v A + u A ) = i-F-. 
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In Eq. (|56|) . F cxt = (F° xt (r, i), Ff^f, t), Fg xt (r, t)) denotes the average external force at point 
(ri,r 2 ,r 3 ) and time t. The vector 



;A(B) s S^K*)*^*) 



(^(f.t),^^*),^^*)) (57) 



is the average velocity of the ensemble, where 

,_ . I 1 if x,- ft) at time t is in a small volume d 3 r around r, . . 

ej(r,t) = < J (58) 

I otherwise. 

Additionally, the osmotic average velocity u A ( B ) appears in Eq. ((56|) . It is connected to the probability 
density p A ( B ) (r, t) Af of finding a particle at point (r, i) in the {A} ({B}) ensemble via 



u 



^Vln(>( B )(r,i)/p ) = ( M f (B) (r,i),4 (B) (r,i),^ (B) (r,t)). (59) 



In Eq. (|59p . po is a constant reference density without physical importance, and v = — — - is the 
friction coefficient of the Brownian motion, which depends on the temperature of the heath-bath and 
Boltzmann's constant. Without loss of generality, v can be written as 

(60) 



2mo 



Since we derive a non-relativistic system, it is possible that for some times, the random and the 
external forces are near infinity. Under these circumstances, a particle may be driven with an infinitely 
large velocity to r = oo at a certain time-step t. However, in physical situations important for the 
non-locality issues discussed later in this article, relativistic effects certainly will occur, and particles 
are not allowed to move faster than the speed of light. Then, an ensemble of quantum particles which 
is emitted at a point [?q, to), for example by a decaying atomic nucleus, can only reach points (r, t) 
within the future cone of (fo,to). Accordingly, we restrict the functions u, v, p(f, t) such that they are 
defined only for these (f, t) in the future cone. 

Before the two ensembles underwent a Brownian motion, it was ensured that their average velocity 
was equal: 

v B = v A . (61) 

It is shown in [15] that this constraint implies u B = u A , too. After the first time step, Eq. ([61]) does 
not hold anymore. Since the {A} particle system gains kinetic energy compared to the {B} ensemble, 
we will get 

\y B {r,t + At)\ < \v A {r,t + At)\. (62) 

We can compute the arithmetic mean 

v(F, t + At) = v A (r,t + At)+^(r,t + At) (63) 

of all particles in the two ensembles. Fortunately, it is always possible to exchange some trajectories 
from the {A} ensemble whose velocity is larger than |v(r, t + At)| with as many trajectories from the 
{B} system whose velocity is smaller than |v(r, t + At)\, until 

v B (r , t + At) = v A (f, t + At) = v(f, t + At) (64) 

holds. 
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The trajectories {A} and {B} are ordered by their velocity and their index j. We may begin the 
exchange procedure by looking at the subset of those trajectories that have the highest velocity in the 
{A} ensemble, and the subset of those trajectories that have the lowest velocity in the {B} system. In 
each of these two subsets, we identify the trajectory which has the highest index j. Then, we change 
the ensembles of the two trajectories that we have found. After this has been done, we check whether 
Eq. (|64|) is fulfilled. In case Eq. (164|) does not hold, we restart the exchange procedure again. The 
trajectory exchange is completed when we arrive at Eq. (|64p . Then, the diffusion process will proceed, 
but for this time with our new division of the trajectories into the two sub-ensembles. 

Every time after the trajectory exchange procedure is completed, we may compute the arithmetic 
mean of the two diffusion equations in ()56|1 . Since Eq. (|64|) then holds, we arrive at 



O 7-1 

— v + v- Vv-u- Vu + z/Au = — v-u- Vu + t/Au= — F oxt . (65) 

at at tuq 

If the Brownian forces were absent, we would have u = 0, with Eq. (|65[) becoming the equation of 
motion of a classical particle. 

For the description of quantum systems, we will assume the Brownian forces to be present. We 
may simplify Eq. (|65|) with the definition of a function 



^(r,t) = ±v^(r^)e^^. (66) 

This function is connected to the ensemble average velocity v via v = — V(f(r, t). (We note that it 
is only possible to formulate the latter equation because one can show v to be curl free.) After some 
computations, Fritsche and Haugk are able to cast Eq. (|6"5"|) directly in form of the time dependent 
Schrodinger equation: 

th^^=U(r,t)^(r,t). (67) 

In Eq. (|67[) . H(f, t) denotes the Hamiltonian operator 

H(r,t)^-^+V cxt (f,t), (68) 

ZrriQ 

with V oxt (r, t) as the potential of the external force F ext = — VV cxt (f, t). The stochastic process for 
the infinitely large ensemble of single particles leads to a differential equation encoding the ensemble 
average velocities v and u in a wave-function. This implies for the interpretation of a quantum me- 
chanical state \4>) that all the information \ip) carries is about the evolution of average velocities from 
a statistical ensemble. It does reveal nothing about individual particles. One should mention that the 
quantum mechanical operators, their eigenvalues and uncertainty relations can be derived from the 
given formalism, as is fully documented in [15]. 

4.2 Quantum states in superpositions 

Before we can discuss the locality properties of the Fritsche-Haugk theory and the physics behind 
entanglement in many-particle systems, we turn our attention to superpositions of two states in a 
single-particle system. The Schrodinger equation is linear. Therefore, the sum of two solutions ip^(r,t) 
and ipi{f,t), i. e. 

*(r,t) = a^ T (r,i) + 6Vl(r,t) (69) 

with coefficients a, b <E C, is also a solution of Eq. (f67|) . With the state in Eq. (|69|) . one would get 
mixing terms in p(r, t): 

p(r, t) = a 2 tf(r, t) + b 2 ^(r, t) + a*b^*(r, t)^ (r, t) + b* a^ (r, t)^(r, t). (70) 
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(72) 



In turn, the Eqs. (|59| would become quite complicated. Accordingly, we need to make a basis change 
in order to get simpler equations. Since both states i/ji(f,t) and ipi(r,t) are solutions of the same 
associated Schrodinger equation, we can write the latter without loss of its physical information in the 
form given below: 

= ih— (6^t(f,t) + oV»;(f,t)) = W,t) (ail>i(r,t) + b^ i (r,t)) (71) 

ifl^(a^(f,t)) = H(r,i)(a^ T (f,£)) 

+ + 

ih^iW^t)) = H(f,i)(6Vi(r,t)). 

Using basis vectors in Dirac notation defined through \ipf) = I _ ) = |+) and \ipi) = I _. 
we get the following coefficients in coordinate representation: 

Cl (r,t) = (^ |tf(f,t)> = ai/>i(f,t) and Ci (f,i) = (^ |*(f,t)) = b^(r,t). (73) 

Hence, the overall state vector \&(r, t)) in Eq. (|69|) becomes 

|*(f,t))=c T ^T)+c i |Vl>, (74) 

and the probability density can be written without mixing terms: 

p(r,t) = {*{f,t)Mr,t)) 

= a 2 ^ (r, t)il>, (F, t) + b 2 ^ (F, *M (F t) 

^^(r^ + fe 2 /^). (75) 

For a 2 and 6 2 , it follows from Jp-|-(4.)(rL,f) d 3 r = 1 and J p(f, t)d 3 f = 1 that a 2 + 6 2 = 1. For a given 
particle in the entire ensemble which is represented by the state in Eq. I|74|l . we have a probability of 
a 2 that the said particle belongs to the sub-ensemble described by ifo (r , t) . Similarly, the probability 
to find that particle in the sub-ensemble represented by ipi(f,t) is b 2 . 

Using the methods from section 4.1, we can derive the two Schrodinger equations for the j and J, 
states in Eq. ([72]) . We start with four Langevin equations 



*B (,\ i m 0,^B r,\ _ T?cxt/,\ , TrBrown/ 



mo^ t (*) + -f^(t) = Yf\t) + F%°™(t), 

m $h(t) ^x£(t) = Ff(t) + F?f°™(i), 

,1 (76) 

mo^(t) + -^x* (*) = Ff(i) + Fff wn (*), 

m ^|(*) ~ ^^(*) = flfW + F^°™(t), 

for the hidden trajectory variables that correspond to the single particles in the sub-ensembles {A |}, 
{B "["}, {A j} and {B J,}. After the same steps that we have already done in section 4.1, we are led to 
the following differential equations: 

-pcxt 

too 

T?CXt 

~ m ° (77) 

T^CXt 

TOO 

T^CXt 

TO 
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dt v T 


Uf ) + (vf 


+ uf) 


• V (vf - uf ) - z/A (vf - uf 




ttf ) + (vf 


-*f) 


V (vf + uf ) + i/A (vf + uf 


-(vf- 


uf ) + (vf 


+ uf) 


• V (vf - uf ) - uA (vf - uf 


9 /-A 


uf ) + (vf 


-uf) 


V(vf + uf)+,A(vf + uf 



The vectors vWn,iIi-(n ar e the ensemble average velocities of the | (J.) states. In this notation, the 
four equations in (|77[) would define two separate ensembles. Every particle would belong to either the 
{"["} or the {1} state for all the time. Nevertheless, for a given single particle of a superposed state, 
there has to be a probability greater than that this particle can be found in each one of these states. 
Fortunately, the two ensembles are no more separate if we demand that some of the fast parti- 
cles in the {A "["} ({A J,}) ensemble are exchanged with the same number of slow particles from the 
{B 1} ({B J}) system. This has to be done with as many particles until 

vf(r,i)=vf(r,i), u*(r,t) =uf (r,t), vf(f,t)=v$(r,t), u B (?, t) = uf (f,t) (78) 

holds. If we add the {A |} and {B J.}, as well as the {A J.} and the {B "f } ensembles, we arrive at two 
equations: 



d_ 

dt 



v? • Vvf - u? • Vu? 



+ vf-Wf-uf-Vu? 



£(*+*?) 



vf • Vvf - uf • Vuf 



vf • Vvf - uf • Vu? 



i/A(u? 



-A(uf 



»>o 



/pext _|_ j?cxt\ 



4) = - ( 

1 mo V 



riext i T?cxt 



(79) 



They can be summed together, yielding: 
d 



-(v B +vf + vf+vf 



vf • Vv* - uf • Vu* 






vf Vvf - uf • Vuf 



vf • Vvf - uf • Vu i 



/A(u| 



uf + uf 



uf = — IF' 



TiCXt | T?CXt 



(80) 



Both equations in (f79"|) contain {A} and {B}, as well as {|} and {J,} particle systems. These equations 
describe the behaviour of equations of motion under the same random and external forces on the 
average. Hence, we have: 



,-?b 



U T = U t = Uj, 



l T 



l T 



(81) 



Using Eq. (|5T)) . we get 2vj = 
Now, Eq. ([50)) becomes: 



*? + vf 



and 2vi = 



-vf , as well as 2u% = u? + uf and 2u = u? 



9 ,_ /_ _, 

— (vf + V|) + I v| • Vvf — U| • Vu| 



+ (v x • Vv 4 - u x • VuJ + ^A (u T + u x ) - — - (F cxt + F cxt ) 



(82) 



This is a sum of two separate differential equations for the ensemble averages of the | and the J, state. 
Its physical information may be written as 



— v T - ut • Vut + ^Au T = F ext , 

at too 

4 vj -3j -Vu, +J/AU, =— F oxt , 



(83) 



r/f 



too 



from which the two Schrodinger equations ([?!?)) can be derived. The two Schrodinger equations in 
(|?2")) can then be added together and one finally arrives at Eq. (fTTj) . The two conditions of Eq. ([75)) 



1G 



and Eq. I|81|) imply that the average velocities vj = v^ = v and uj = u*j = u are equal. Hence, in 
the interpretation of stochastic mechanics, Eq. (|71[) describes two oscillating ensembles with identical 
dynamics on the average. 

The procedure above also works with an arbitrary number of states \ipi(r, £)), \ip2(F,t)), ■ . ., 
\i/>fif(r, t)) that are summed up to a general superposition with weighting factors aj G C: 

M 
\^(P,t))=J2^Mr,t)). (84) 

i 

The weighting factors a, define the probability of a given particle to be found in the i-th state \ipi(r, t)) 
of the superposition. Therefore, the ensembles in the corresponding J\f diffusion equations must contain 
different numbers of particles. More precisely, each state \ipi(r, £)) has to be computed with an ensemble 
of Saf trajectories j\,J2,J2, ■ ■ ■ iJsa 2 - This ensures that there is a chance of af to find a particle in the 
system associated with \ipi(r,t)). 

Furthermore, the exchange procedure must be defined differently. Two trajectories with indices 
j and f are now exchanged between two ensembles denoted by {Bi} and {Afc}, where i,k £ N are 
numbers that range from 1 to TV. The indices i and k define the ensembles corresponding to the states 
\ipi(r,t)) and \ipk(r, i)) between which the trajectories are interchanged. We note that i and k have to 
be chosen randomly, but with equal probability for each state. The different probabilities of finding a 
given particle in the separate states \ipi(f, £)) and \ipk(r 7 t)) only arise due to the unequal size of their 
corresponding ensembles. 

After the necessary interchanges took place, one may add the TV differential equations in one step 
and split them up into separate equations for the 1, 2, 3 . . . , Af systems, analogous to Eq. ([53")) . We 
remark that because of Eq. ([59]) . nothing about the weighting factors af in the probability density 
appears in the ensemble averages v and u. Hence, one has to re-introduce this information after one 
gets from the separated system in Eq. (|72|) to the superposed Eq. ([7Tj) . 



4.3 Locality in the theory of Pritsche and Haugk with single-particle sys- 
tems 

In this section, we turn to the locality properties of the stochastic process described above. We observe 
that in the Langevin equations (|53[) . there is no dynamical coupling between different trajectories 
xf T m(*) and xf^l (t), where j £ f . The force field F° xt (t) only acts on the j-trajectory, and Ff Iown (t) 
is required to be independent for each j. Nevertheless, the time evolution of different trajectories will 
be correlated because of the trajectory exchange mechanism. 

We will now assume that something happens to a large set of trajectories at a point (r',t'). For 
example, the particle speed may be altered there by external forces controlled by the experimenter. In 
case (r ', t') lies in the past cone of another point (r, t), some of the changed trajectories may arrive at 
(f, t). We will denote the changed trajectories that reach (r,t) by xi,:,(t). 

Let x-},,l(t) be another trajectory at (f, t). The exchange procedure decides whether x.. (| (() will 

change its ensemble. This decision depends on the magnitude of the particle's own velocity x.. (| (t) 
and the velocity distribution of all trajectories from the {A} and {B} ensembles that have arrived at 
(f, t), including the trajectories x.AA(i) which were changed at (f',t'). 

For the next time step, dependent on the trajectory exchanges required for Eq. (|64|) to hold, x J,. ( (t) 

might be placed in the {B} ({A}) ensemble. In case x.A, (t) changes its ensemble, this trajectory will, 
at least for the next time-step, be subject to the differential equation of the new ensemble, and we will 
have to rename it as x.},,l(t). 

Since the interchange procedure explicitly depends on the modified trajectories x.L >({), the par- 
ticle trajectory x\,l(t) may be influenced at (f,t) by what has happened at (f',t'). Moreover, for 
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that kind of influence it is not required that x.} ( .l(t) ever went through the point (r',t r ). This effect 
is important in double slit experiments, where a particle's trajectory, after the passage of one slit, 
depends on whether the second slit is open or not, no matter if the particle actually went through the 
second slit. 

However, it is important to note that the trajectories xi,, (() can only be influenced at (r,t) 
through the exchange mechanism if (f',t'), where some intervention was made, lies in the past cone 
of (r, t). In turn, there does not exist any violation of active locality. 

By writing this, one should be aware that we are still discussing a non-relativistic system. The 
external and Brownian forces were given in a non covariant form. Hence, these forces are still able 
to modify the velocity of the particles instantaneously, even though we have restricted the particle's 
maximum speed to lightspeed in section 4.1. This last problem could only be resolved with a covariant 
formulation of the Fritsche-Haugk theory. Fortunately, as we will see in section 4.3, the discussion of 
entanglement does not depend on that problem. The trajectories corresponding to different particle 
states of a many-particle system will turn out to be completely uncoupled, at least if inter-particle 
forces are not present. 

Finally, we emphasise that the exchange procedure should not be seen as a dynamical coupling 
mechanism between different trajectories. Instead, the comparison of a particle's velocity with the 
velocity distributions of the {A} and {B} ensembles is merely a model to account for the idea of kinetic 
energy conservation on the average in a stochastic theory. That there are in fact no dynamical couplings 
between the individual trajectories of the ensemble can be seen experimentally in self-interferometry 
experiments. The typical example is neutron interferometry [32], where "usually, at one time, there 
is only one neutron in the interferometer, if at all because at that time, the next neutron has not yet 
been born and is still contained in the uranium nuclei of the reactor fuel" [33] . 

4.4 Derivation of the many-particle Schrodinger equation 

In this section, we want to investigate whether active locality persists when we consider many-particle 
systems. The N-particle Schrodinger equation can be derived if the three dimensional coordinate space 
of the single-particle problem is expanded into a 3N dimensional space. We start with the Langevin 
equations of the j-th particle trajectory in the ensemble of the n-th particle state, where n € N and 
1 < ?i < N: 



n ^n 



(85) 



In Eq. (|85|) , there is a sum of additional force terms F°*' ( x „ (t) , x ,, n , (t) , t ) corresponding to inter- 
particle forces between the particle states n and n' ^ n. The inter-particle forces are defined to be 
equal for the {A} and {B} ensembles: 

F«? (5f n (t),^,„,(t),i) = F^ (Z%(t),Zf, n ,(t),t) . (86) 

From the two Langevin equations (|85p . Fritsche and Haugk derive the following differential equa- 
tions for the {A} and {B} ensembles, completely analogous to the single-particle case: 



... ,,..„ . , A*NB , ,jNB\ Y7N /,-^NB -NB\ „ A N /,-?NB ,^NB\ _ ^N,ext 



V N (v NB - u NB ) - *,A N (v NB - u NB ) = F f 



d_ 

Of 

9_ ^ NA + -NA) + (.NA _ fl NA) . yN ( ? NA + gNA) + ^ A N (^NA + -NA) = pN,cxt_ 



(87) 
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In contrast to the single-particle case, F N > ext is a 3N dimensional vector of the form 



4- V , 
mo, n=i Z-fT» =1 



F° x j N (^=N,t) , V-N Fg N (f„ =N ,f 



Ely 
n' — 
n ^n 



m ,„=N 



Its n-th component F^ ,oxt consists of the average single-particle force F° xt (r„, t) and a sum of average 
inter-particle forces F° xt (r n , r n >,t). Both forces F® xt (r„, i) and F° xt (r„, r„', i) act on the n-th particle, 
which has rest mass of m 0n . 

The three dimensional ensemble average velocities u and v of the single-particle case are similarly 
replaced by 3N dimensional vectors. They depend on the 3N coordinates r N = (fx, r 2 , . . . , ^n) of the 
many-particle problem in the following way: 

v NA ( B )_(vr ) (n li ),...,vr»(f Nli )) , S NA ( B > ^(4 (B \n,t),...,4 (B) (r N ,t)) . (89) 

The n-th component of u NA( - B ^ and v NA ^ B ' is computed only with the trajectories x ■„ (t) of the n-th 
particle state. This gives rise to the coordinate dependence shown in Eq. (f59"| . 

In order to derive the many-particle Schrodinger equation, an interchange between some trajectories 
corresponding to slow particles from the {B} ensemble and the same number of fast particles from the 
{A} ensemble is necessary after each time-step, until 

v NA (f N , t) = v NB (f N , t)^v N = (vi(rl, t), v 2 (f 2 , t), . . . , v N (r N , t)) (90) 

holds. Here, we have to note an additional subtlety that does only arise in many-particle systems. It 
has to be ensured that exchanges between two trajectories x^„ (t) and x-, n , (t) are confined within 
the same particle state n = n'. An interchange between trajectories of different states n and n' ^ n 
would couple ensembles corresponding to different particle states, even if no inter-particle forces are 
present. This would cause a general impossibility to define single-particle operators, which act on 
only one state n of a many-particle system. Fortunately, due to the specific coordinate dependence 
of the ensemble average velocities v NA ( B ) and u NA ( B ) in Eq. (|89|) . the components of Eq. (j90"|) must 
agree for every state n separately. Therefore, the exchange procedure can well be confined within each 
single-particle state only. 

Every time after the trajectory exchange is completed, one can compute the arithmetic mean of the 
two differential equations in (|87| . Analogous to the single-particle case, we get the diffusion equation 

-v N - u N ■ V N u N + ^A N u N = F N ' cxt . (91) 

at 

From Eq. (|9ip . the many-particle Schrodinger equation 

5* (r N t) 
ift ^ ' J = ft (?*,*)* (?".*) ( 92 ) 



can be derived. In Eq. (|92p . the Hamiltonian operator is given by 

-h 2 A] 

2m n 



N ^2 AN 



H (<*.*) ^£l^+ V «* (**'*)■ ( 93 ) 



n=l 

Its external potential 

N N N 

V oxt (r N , t)=Y, V cxt (f„, t) + - Y.Y. Vext(rk, f n >, t) (94) 

n=l n=l n '=i 

n ^n 
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is connected to the n-th component of F N,cxt 



and the velocity distribution of 



via 
F^ = -^-V^ t (f N ,t). (J):,] 

4.5 Locality in the theory of Pritsche and Haugk with many-particle sys- 
tems 

In this section, we assume all inter-particle terms to be absent. We observe that with the absence of 
inter-particle forces like F|jf (x „ (i),x/ n , (£),£] , there are no terms in Eq. ([55]) which interconnect 

different trajectories x „ (t) and x-,„> (i) from separate particle states n' ^ n. A similar statement 
is also true for different components of the ensemble average velocities v^ and u^, because they are 
calculated with the x -J ; ' (t) trajectories of each single-particle state n separately. 

The decision whether a trajectory x ■ (t) changes to x -^ (£) during the interchange procedure 
only relies on the magnitude of the (J, n)-th particle's velocity x ■ (t) 

all particles that belong to the state n and are in the past cone of x n (t). Therefore, the trajectories 

x-J '(t) are never coupled in any way to the trajectories x -,„,(£) of other particle states n! ^ n. As a 
result, we conclude that active locality holds in the hidden variable theory of Fritsche and Haugk. 

4.6 Entangled states in the model of Fritsche and Haugk 

That there does not exist any dynamical coupling between the different particle states which a many- 
particle system is composed of doesn't imply that correlations between these states are forbidden. In 
section 6, we will introduce further conditions on the Gaussian distributed Brownian forces F^ own (£) 
in a completely local framework. Those conditions will be such that under certain circumstances, the 
trajectories x„ (t) change their ensemble from {A} ({B}) to {B} ({A}) simultaneously with the 

trajectories x /„,(£) from another state n' ^ n. Then, correlations between ensembles corresponding 
to different particle states will emerge. The resulting many-particle wave-functions are of the form 

\*(r u r 2 ,t)) = -= \MriA) <8> \M^, *)} + ^= \Mn,t)) <8> \MM) » ( 96 ) 

and are commonly called entangled. They were observed for the first time in massive fcrmionic systems 
with a pair of entangled protons in 1976 [34]. 

The state vector ^(fi, r*2, t)) in Eq. (|M|) is an element of the product Hilbert space TL = Hi (8> 7^2- 
It is a superposition of the two product states 

\Xi(n,f 2 ,t)) = \Mn,t))^\Mr2,t)) and \x2(T f i,T f 2 ,t)) = \il>2(* f ut))®\Mi f 2,t)), (97) 

where \ipi(?i,t)) , \ip2(n,t)) e Hi and \ijji(r2,t)) , |V>2(j2, t)) G H 2 - 

The entangled state solves a many-particle Schrodinger equation with an Hamiltonian operator 

2 _t2AN 2 

H(r 1; f 2 , t) = Y, ~ " + E Voxt(r„, t), (98) 

n=l 0n n=l 

which has no inter-particle forces and whose external potentials agree at each time t: 

V ex t(f ? i,t) = V ext (r 2 ,t). (99) 
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The Schrodinger equation associated with ^(ri, rb, t)) can be written, analogous to Eq. (f72|) . as: 

ih— |*(ri,r 2 ,t)) =H(ri,r 2 ,t)|*(ri,f 2 ,t)) (100) 

' lh m (72 IXi(n,r 2! t))) = R(r u r 2 ,t) (^ |xi(^,^,*)>] 

+ + (101) 

^I(tI IXa(n,^, *))) = H(r 1 ,r 2 ,t)(^| X2 (r 1 ,f 2 ,i))). 

This system consists of two particle ensembles r\ and r 2 . The particles from fi can be in different 
states IV'i(fi)) or 1^2(^1)), whereas the particles from r 2 can be in the states 1^1(^2)) or |^ 2 (r 2 )). 
In order to describe the dynamics of the r\ particles, the trajectories in the {Al} ({Bl}) ensemble 
of IV'iCn.)) must be exchanged with the trajectories of the {B2} ({A2}) ensemble of \ip2(?\)) 1 and a 
similar procedure has to be done with the r 2 particles. In the meantime, one has to ensure that each 
j'-th pair of particle trajectories in the r% and r 2 systems changes its state simultaneously. We will 
show in the context with Bohm's EPR experiment how this can be achieved although the trajectories 
of the two different ensembles are disconnected from each other. For now, we only remark that without 
inter-particle-forces, it is in fact a requirement for the particle trajectories x -^ (t) and x-,^ (t) not to 
be interconnected directly in any way. This may also be seen with the following example by Nelson [8] : 
If one induces an external measurement process on the n-th particle state of a multi-particle system, 
one may get information about the hidden variable x ■ (t) . The measurement device may change 

x n (t) in a way that correlations with all the other particle trajectories x-, n , (t), where n' 7^ n are 
lost. Since the different particle systems are uncoupled, an interaction happening on the ensemble of 
the n-th particle state can in no way influence the properties of the particles in the other states nl . 
Let H = Hi tg) H.2 be the product Hilbert space above and define an Hamiltonian operator 

H = H x <g) 1 2 + li <g> H 2 , (102) 

where Hi € Hi and H 2 s 7i 2 . The time evolution for any observable Ai e Hi in the Heisenberg 
picture 

A ft (t) = e^Axe- 4 *" = e itftl Aie"^ 1 ® 1 2 (103) 

is completely independent of H2, even if |^(fi, r-z, t)) € H is described by the entangled state of Eq. 
fll|. As Nelson put it: 

"No matter how systems may be correlated, there is no way of telling, what [external] forces 
the second system is being subjected to, as long as they are not dynamically coupled." 

In this context, we want to mention that Nelson's original theory could not reproduce the example 
above [8] because he identified v and u not as average velocities of an ensemble but assumed them 
to be the trajectories of individual particles. With this assertion, neither the splitting of the many- 
particle Schrodinger equation (jTUTJ) , nor any sort of exchange procedure between the particles in the 
{Al} ({Bl}) and {B2} ({A2}) system was possible. Accordingly, active locality gets lost in Nelson's 
original formulation of stochastic mechanics, similarly as it gets lost in Bohmian mechanics. There, 
(see [62] for a general introduction) it is well known that the derivation of the many-particle Schrodinger 
equation is only possible provided that a highly non-trivial dependency between the Bohmian "tra- 
jectories" of different particle states is introduced, which is at odds with the theory of relativity and 
without any physical foundation. Interestingly, Nelson himself has shown recently [35] that v and u 
cannot correspond to properties of individual particles since this would lead to measurable differences 
between his original theory and ordinary quantum mechanics. A similar observation was made for 
Bohmian "trajectories" earlier [36]. 
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Figure 3: A particle with spin | rotates in the left direction with respect to the z axis. 

5 Spin in a Fritsche-Haugk like hidden variable theory 

5.1 Preliminary considerations 

Nelson's derivation of the Schrodinger equation was quickly extended to include particle spin, but 
without proper physical interpretation [37]. Over the years, more physically motivated derivations of 
the Pauli equation were developed. One by Faris [38,39], and another one by Fritsche and Haugk 
[40]. Unfortunately, these frameworks did either violate active locality or they did not contain any 
mechanism that would change a particle's spin state after its emission. Thereby, none of these theories 
led to a violation of deterministic passive locality only. In this section, we propose a theory of spinning 
particles providing such a mechanism. 

Beginning with an example of [40], we try to interpret the spin of a particle as rotation around 
a chosen axis z (see Fig. [3]). We let the particle rotate around this axis under energy fluctuations 
AEAt i=s h/2. For the rotational energy, we have the expression E = ^-cOq from classical mechanics. 
In this equation, O = m^r 1 denotes the inertial tensor with itiq as the particle's mass and r as the 
distance of the particle from the chosen axis z. The particle's angular velocity is given by ojq = ^ > 
where T is the time needed for the particle to perform a complete rotation. 

If we identify At with T and AE with E, we may write 



a r.A 1 27r27r 

AEAt = —m$r ■ r T = m^vrix — ns rs a/2. 



Ineq. (HID, 



movr 



27D- 

T 

h 

' 2^ 



is the magnitude of the particle's velocity and 



(104) 



(105) 



defines the magnitude of its angular momentum. Apparently, s is independent of the particle's mass 
mo and the radius r of the rotation. The arguments above were only approximations. For later, it will 
be useful if we set the magnitude of the spin angular momentum to s = h/2 on the average. 
From the time dependent Schrodinger equation, one can derive the theorem of Ehrenfest: 



d 



^ v = F 



(106) 



In Eq. (|106p . Mis the expectation value of the velocity operator and /F ext \ is the expectation value 

of the external force. Ehrcnfcst's theorem implies that the angular momentum operator L = r x m^v 
may not vanish in the case of ft — > 0. On the other hand, the quantum mechanical spin should always 
go to zero in the classical limit by Eq. (|105|) . Hence, the particle's spin must be completely independent 

of L and we must expand our description. 

5.2 Stern-Gerlach effect 

If particles of spin 1/2 are not specially prepared, they occur in a superposition of | and J, states. 
This superposition can be derived from a similar computation as in section 4.2. Here, however, the 
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trajectories x.},.l(t) of Eq. (|7S|) describe particles in rotational motion. Their rotation can be in the 
left or in the right direction with respect to the chosen axis z. Particles rotating in the left direction 
are denoted by j, whereas particles rotating in the right direction are identified by J.. After some 
calculations, we eventually arrive at the four differential equations of Eq. (177|) . 

As explained in section 4.2, a given particle in an ensemble described by a superposed quantum state 
must have a non-vanishing probability to be found in each one of the component states. Accordingly, 
we have to exchange trajectories between the {A "["} ({A J,}) and {B J.} ({B ]"}) ensembles before we 
get to the differential equations (|79|) . A nice drawing of this process is given by R. Penrose [41] for 
Dirac particles. It is this procedure, whereby the possibility for a violation of deterministic passive 
locality does arise. 

At the end of the calculation in section 4.2, we were able to describe the behaviour of the superposed 
system by a probability density 

p(f,t) = a 2 pi(r,t) +b 2 Pi (f,t), (107) 

and a quantum state 

|*(f,t)> = ( l^p^ ) =a^{r,t)\+)+b^{f,t)\-). (108) 

Both are in form of a sum of entities from the two rotating particle ensembles. 

If an external magnetic field Bg(r, t) is applied in the direction of z, the energy of i/j^ (r , t) increases 
with a magnitude of fXsBg{f,t) \a\ , whereas the energy of ip\(r,t) decreases by /iBBg(r,t) \b\ . The 
constant fiB is called Bohr magneton 

*■- 55- (109) 

It emerges since the relativistic Dirac equation is, in contrast to the Pauli equation to be derived 
here, linear in all its derivatives. The necessity of this linearity, which is absent in the non-relativistic 
theory, changes the interaction of the angular momentum with magnetic fields. If one does linearise 
the Schrodinger equation in a similar manner as it is typically done for Dirac's equation, one gets the 
same value for \iq [42]. 

The changes of the quantum mechanical energy densities, which are caused by the interaction of 
the spin rotation with the magnetic field Bg(r, t) , are given by the expressions 



and 



^{r,t) =^i B Bg(r,t)\a\ 2 ^Ur,t)^{f,t) (110) 



/*;(?, t) = -MB^(r,i)|&| 2 ^(r,^j.(r,£). (Ill) 



With their sum 

p{f,t)=^{r,t)+ H {r,t), (112) 

the total change /j,(r, t) for both spinors can be cast as 

H(f,t) = *+(r,t)^B^(f,t), (113) 

where we have used the definition 

B*=B*(n*)( J ^ )• (114) 

The different coupling of the i>ui)(r,t) states to magnetic fields gives rise to the Stern-Gerlach effect 
[43] . If atoms from an incident beam with equally distributed magnetic energy density enter a strong 
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magnet, the two spin systems are subjected to different forces F°5 t (i) 7^ F°j ct (t). Accordingly, they get 
deflected in opposite directions and we are left with disconnected trajectory ensembles corresponding 
to the different spin states. Due to active locality, the exchange procedure can only happen between 
trajectories x A , I (t) that are located at the same point (r, t) in space-time. Hence, after the separation 
by the Stern-Gerlach magnet, the exchange procedure can only take place between trajectories of the 
same spin ensemble. This is observed in double Stern-Gerlach experiments [44]. If one places a second 
Stcrn-Gerlach magnet behind the first and let it have the same directional orientation, the parts of 
the beam that were previously separated, do not split up again. In this way, particle ensembles in a 
beam can be prepared to contain one spin state only. 

The repeated Stern-Gerlach experiment demonstrates both the locality and the non-markovian 
property of the stochastic process involved. The stochastic differential equation that a particle is 
subject to depends, due to the coupling of different trajectories in the exchange procedure, on what 
has happened to the entire ensemble before. Therefore, the stochastic process that acts on a single- 
particle has a kind of memory about the past interactions between the magnets and the ensemble in 
which the individual particle is currently placed in. Nevertheless, the stochastic process is a local one. 
After the two spin systems are spatially separated, they do not influence each other anymore. 

5.3 Pauli equation 

From here on, we can closely follow [40] to derive the remaining parts of quantum mechanics for a 
particle with spin. In [15], the Schrodinger equation under a magnetic field with its Hamiltonian 

(-iKV - eX. (r, tj) 

W,t) = ± J —+e${r,t) + V cxt (r,t) (115) 

Zrrio 

was derived from a Brownian motion process. If the two rotating systems undergo this process, their 
overall energy density maybe written as superposition 

w(r,t) = tf(?,t)1[(f,t)M?,t) +1>t(f f ,t)R(f,t)il> l (f,t) = *+(r,i)HvI,(r,i), (116) 

where we have used the definition 

H = H(r*,t)( I J V (117) 

The Hamiltonian should include the energy from the interaction of the spin rotation with Bg(f, t) . 
Hence, /J,(r,t) must be added to w(r,t), and we get H + /UbB^ for the Hamiltonian operator. We 
observe that the effect of ^B| on each spin state is that of a potential. Since the time dependent 
Schrodinger equation was derived under an arbitrary external potential, we can write a first equation 
for a spinning particle: 

iH ^dt^~ = (** + ^ bBz ) * (r *' t] - (U8) 

If we want to measure the particle's spin rotation around an axis z' that is different from the axis 
defined by the z component B z (r , t) of the magnetic field, we must change the coordinate system of 
B(r, t). This can be done via linear transformations B'(r, t) = QB(r, t), where Q defines a 3x3 rotation 
matrix which transforms (B x (r,t),B y (f,t),B z (f,t)) into B'(r,t) = (B X '(f,t),B y i(r,t),B z >(f,t)). It is 
known from classical mechanics (see e.g. [45]) that by using Euler angles ijj, (f> and 9 as shown in Fig. 
01 any 3x3 rotation matrix can be converted to the following 2x2 form: 

/ e^ +< « cos (6>/2) ie*W-0sin(0/2) \ , , 

^-\ie-^-*Um{d/2) e-*<*+*>cos(0/2) )' [ ' 
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Figure 4: Two coordinate systems with unit vectors (S,y,z) and (x' ,y',z') related by Euler angles 
r/>, and 6*. The angle tp = <fi — ir/2 is also depicted 



In order to apply the 2x2 matrix from Eq. (|119|) on B(f,t), we must rewrite the magnetic field: 

b= i „ ,-/;:;',".;: , _, ,. ""'•';, , ir// i . (120) 



B z {r,t) B x {f,t)~iB y (r,t) 

B x (?,t)+iB y (r,t) -B z {r,t) 



Then, one can compute B' via 
B' = QBQ+. 



(121) 



We note that a magnetic field B(r, t) = (0, 0, B z (f, t)) has a matrix form which is equal to B^ in Eq. 
I|114|L If we use the Pauli matrices a$, ag and ag, we can simplify equation (|120|) considerably. At 
first, we observe that 



B = B x (r, t)a s + B y (r, t)a$ + B z (r, t)a s . 
Writing the Pauli matrices in form of a vector a = (ag, ag, ag), we arrive at 
B = B(r, t)a. 



(122) 



(123) 



It has to be derived, how ^(r, t) will transform under a rotation of the coordinate system. At each 
point (r, t), the change of the particle's energy density, which is caused by the interaction of the spin 
with the magnetic field, should be invariant under a rotation: 



Making the assumption that 

*'(f,t) = Q*(r,i), 
we get the following expression for //(r, £): 

M '(f, t) = ^ B * + (f, £)Q+B'Q*(r, t). 
Since Q in Eq. (|119[) is unitary, the following identity holds, provided we use Eq. (|121[) : 

Q+B'Q = Q + QBQ+Q = B. (127) 



(124) 
(125) 
(126) 



The equations (|127[) and (|126p imply that the condition of Eq. (|124|) is fulfilled. This justifies our 
assumption from Eq. (|125p . 

Using Eq. (|126p , Eq. (|127p and equation (|123p . we can write the overall energy density of a particle 
ensemble under a general magnetic field: 

(128) 



w (r i t) = #+(f, t) ( H + HbB{P, t)a *(f , t). 
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Again, the energy from the interaction of the spin with the magnetic field has the form of a potential. 
Hence, we finally arrive at the time dependent Pauli equation for an arbitrary magnetic field: 

iH ^dt^ = {^ + ^B(f,t)a) *(r,i). (129) 

In order to describe the direction of the spin rotation around the axis z' , only the Euler angles 9 
and <f> are necessary, and one can select ip — —ir/2 without loss of generality. The projection of the 
rotated axis z' in the z/y plane of the non-rotated coordinate system encloses an angle ip with the x 
axis. This angle is related to the Euler angle <p by <j> = <p + 7r/2 (see Fig. [4]). Inserting ip + ir/2 and 
the choice of ip = — n/2 in the matrix of Eq. (|119[) usually simplifies the calculations. 

A common situation is that the magnetic field B(r, t) of a Stern-Gerlach apparatus lies in a coordi- 
nate system which is rotated with respect to the orientation of an incoming spinor \t'(r, t). Especially 
in EPR experiments, one has to compute *f>(r,t) = Q + \t'(r, t). Applying the adjoint of Q with the 
aforementioned simplifications on a ^'(r, t) = ipi(r, t) |+) spinor yields 

y(r,t) = ^(f,t) (e^/ 2 cos (0/2) |+) +e 1 ^ 2 sin (6/2) |-)) , (130) 

and we get a similar result for a ^f'(r, t) — ip'Af, t) |— ) state: 

*(r, t) = i>{ (r, t) (-e-^' 2 sin (6/2) |+) + e"»/ 2 cos (6/2) |-)) . (131) 

5.4 Spin operator 

Still following the lines of [40], we will derive the so-called spin-operator. Both rotating particle 
ensembles have a spin angular momentum of ±H/2 on the average. Therefore, the overall expectation 
value of the spin in the z direction can be written as weighted sum 

(Sg) = ^J (\a\ 2 ^(r,t)^(r,t) - |6| 2 ^(r,t)^(r,t)) d 3 r. (132) 

We can simplify this expression with the definition of a matrix 

fc-KJ-i)-^ (133) 

and obtain 

(£*>= [iP+(r,t)S^(r,t)d 3 r. (134) 



For the spin with respect to the direction of an unit vector z' = jl from an arbitrary rotated coordinate 
system, we get, analogous to Eq. (|123|) . the relation 

S?' = ^pa. (135) 

Equation (|I35[) is the so-called spin operator for the z' direction. It is convenient to write the axis 
z' in spherical coordinates as z' = (cos(psm6,smipsmd, cos 6). We observe that jl depends on the 
same angles that are used for the spinors of Eq. (|I30[) and Eq. (|I3f[) . Finally, the average value (Sgi) 
becomes 

&>) = [i> + (7,t)§z^(r,t)d 3 r =ft(Sg) , (136) 
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where 

%) = /V(r,i)S^(r,£)d 3 r. (137) 



Having derived the Pauli equation, the correct transformation properties of the spinors in rotated 
magnetic fields, and finally the Spin operator for arbitrary magnetic fields, we should look again at the 
Stern-Gerlach effect. Usually, it is said in textbooks that a Stern-Gerlach magnet would simply be a 
measurement device for eigenvalues of the Spin-Operator. However, a closer look shows that this is an 
oversimplification, that would only work for an infinitesimal small Stern-Gerlach magnet. In the real 
world, the magnetic field of a Stern-Gerlach magnet is intricately distributed over a broad part of the 
space. In order to take this into account, we have to solve the Pauli equation with a rather complicated 
magnetic field. One then finds that the results may differ from the naive eigenvalue computations. 
Such effects are especially important in experimental setups where the particle beam successively goes 
through two Stern-Gerlach magnets. The reader might consult the articles [46] and [47] for further 
information on this issue. 

6 Entangled systems with spin 

6.1 Transformation of the entangled state into equations of motion 

As in the case of the many-particle Schrodinger equation, we can extend the coordinate space to derive 
the N-particle Pauli equation. From now on, we will assume that there are no forces which interconnect 
the different particle states. 

We start our discussion with a state 

^(r ll r 2 ,t) = ^=(Mri,t)\-)^Mr2,t)\+)-Mn,t)\+)^Mr2,t)\-)) 

= ^Xi(n,r 2 ,t) + -^ X 2(ri,r 2 ,t). (138) 

It shall describe two ensembles of entangled particles with the same rest mass mo and opposite spin. 
The one-particle spin states in the first line of Eq. Q138|) are all eigenstates with respect to a chosen 
axis z and in the second line, we have used the functions 

Xi(ri,f 2 ,t) = Mn,t) |-) ® ip 2 (r 2 ,t) |+) and X2(n,f 2 ,t) = ifafat) \ +) <g> fa (Fa , t) |-) . (139) 

The state in Eq. (|138[) is the solution of a two-particle Pauli equation without external electromagnetic 
fields. This leads to a Hamiltonian operator 

/ 2 -fr 2 A N 2 \ / 1 n \ 

H(*,ft,t)= (E-2^ "+£ W*,*)J ( 1 ) ' (140) 

whose external potentials for the two particle systems agree at each time t: 

V ext (fl,t)=Vext(f : '2,t). (141) 

Similar as in Eq. (jlOlj) . the Pauli equation for the entangled state in Eq. (I138J) can be written as: 

ih— ^(fi,f 2) t) = H(r 1 ,r 2 ,t)V(r 1 ,r 2 ,t) (142) 



; /'4r f^Xi(fi,f 2 ,i)j = H(fi,f 2 ,t) (^Xi(n,r 2 ,t)J 



"at 



^V 



(73X2(^1, r 2 ,t)j = H(fi,F 2 ,<) {^X2(ri,r 2 ,t)J . 



(143) 
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The probability density of the particle ensembles described by Eq. (|143|) is given by 

p{ri,r 2 ,t) = -pi{fi,f 2 ,t) + -p 2 (ri,r 2 ,t), (144) 



where 



Pi(ri,r 2 ,t) = xt(ri,r 2 ,t)xi(ri,r 2 ,t) 

^^l(ri,t)4> 1 (ri,t)h\-)®^ 2 (r2,t)Mr2,t)(+\ + ) 

= p n (ri,t)®P2i(r 2> t) (145) 



and 



P2(ri,r 2 ,t) = xt(ri,r 2 ,t)x 2 (ri,r 2 ,t) 

= ^(ri,t)^2(ri,t)(+| + )®Vi(^2,i)^i(r 2 ,t)(-|-) 

= p2t(fi,t)^p n {f 2 ,t). (146) 

The osmotic velocity is connected to the Probability density of many-particle states via 

u N (,- N ,0--^V N ln(p(r N ; 0/Po). (147) 



Hence, after the tensor product in pi(fi,f 2 ,t) is converted by the logarithm laws into a sum, we get 



for ut^ri,^,*): 



-'"-" ; ■'"■' ''"' I - I ",■■>>' \ ( 14 g) 



u^(r 1 ,r 2 ,t) = -^-V N (ln((p u (f 1 ,t))/p )+ln((p 2T (r 2 ,i))/p )) 
/ r * tf]n((pu(rM))/po) \ / tt|(fi,t) 

u T (ri , i) 
, (,; id ui unc lunu u x2 v > <-; = I 

The definition of the average velocity v (r , t) for a many-particle ensemble is 



Similarly, uS, in the Pauli equation for x 2 (r*i,r 2 ,t) is of the form £ (r N ,t) b. , 

V ui(r2;t) 



v N (r N ,t)^— V N ^(r N ,t), (149) 

where y> (r N , t) is the phase of the many-particle wave- function 

* (f N ,t) = ±^Jp{r N ,ty v ( r ~ N ' t ). (150) 

Writing xi(r*i,r 2 ,t) in the form of Eq. (|150[) . i. e. 



Xi(ri,r 2 ,t) = sj pn{n,tym^ ® ^j^ e ^( ft ''), (151) 

we can identify the overall phase of Xi(^ii ^2, t) as 

^i(ri,f 2 ,t) = ipii(ri,t)+ip2i(r 2 ,t). (152) 

Using the Eqs. (I149|) and ||152|) . we arrive at 



: -\ ,-■ ... , , _ S yvN,„ /-. -. .s_( n(ri,i) 
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(154) 



and we get v^ 2 (r*i, fi, t) = I _ ,., ' » I in completely the same way. 

The expressions for u^ , u?? 2 , vrfj , v^ 2 imply that the Pauli equation (|143[) with the entangled state 
in Eq. (|138[) can be decomposed into 2 two component differential equations, where F 2 xt (r 2 ,i) = 
Ff^fi,*): 

| v x (rl , t) - V h (n , t) • Vu x (n , t) + z/Aui (rx , t) = ^Ff * (* , t) 
|v T (f 2 ,i) - u T (r 2 ,£) ■ Vu T (r 2 ,£) + j/Au T (r 2 ,i) = -i-F c 2 xt (r 2 , £) 

f t v T (fx , t) - S T (n , t) ■ Vu T (n , t) + i/Aift (n , t) - ^Ff * (f x , t) 

^vx(r 2 ,t) - u x (r 2 ,i) • ^(fa,*) + i/Aui (f 2 ,t) = i^(5,t) 

The physical information of Eq. (|154|) is equivalent to what is described by Eq. (|143[) . However, 
Eq. (|154p contains the ensemble averages u and v of each particle system directly. All equations in 
(|154[) can be similarly derived as the two equations in (|55|) . but the derivation must be done for the 
r\ and r 2 particle ensembles separately. 

6.2 How the correlations of entangled states emerge in the local hidden 
variable theory 

Having converted the entangled state into equations of motion that can be derived from the previously 
discussed hidden variable model, we will investigate in detail how the correlations of entangled systems 
emerge from the viewpoint of these hidden variables. 

The source of the EPR experiment produces particles of a defined spin direction with respect to 
the axis z. However, the source cannot prepare the spin states in separate beams like a Stern-Gerlach 
magnet because it does not divide the two spin ensembles spatially. Hence, the spin states of the 
particles emerging from the source are in a superposition of an ] and a J, state. According to section 
5.2, the individual particles in such a system have to undergo an exchange procedure between the 
{A |} and {B J.} or the {B ]} and {B J,} ensembles. The trajectories x.i,, L(f) of the f\ particles have 

to switch between their vf(fi,t) and v?(fi,t) system, as well as between their v£(ri,t) and vf (fi,t) 
ensemble in Eq. (|154p during their flight. An analogous procedure has to happen with the trajectories 
xJv,> 2 (£) of the r*2 ensembles. 

We observe that it is possible, to threat the r\ and r 2 particle systems completely disconnected 
because they are uncoupled in Eq. (|154p . These ensembles may even be spatially separated, in a 
way that there is no contact between them. A local measurement device that acts on one particle 
ensemble, say r\, might change the potential V ex t(ri,i) in Eq. (|140p . Along with this goes a change 
of the external force F° xt (r*i , t) . If the action of the device was such that Ff xt (r*i , t) ^ F 2 xt (r 2 , t) at a 
certain time t, the equations governing the f± particle system would be different from the equations of 
the r 2 system. As a result, correlations in Eq. (|154[) between both ensembles would be lost. Obviously, 
the dynamics of the r 2 ensemble is not be changed by something that happens with the r% ensemble. 

For the discussion of the EPR experiment, it is important to note that there will still be correlations 
between the two particle systems if the experimenters at both measurement stations 1 and 2 choose 
the settings of their devices in a way that Ff xt (fi,£) = F 2 xt (r2,i) for all times t. This is exactly the 
situation when the two separated Stern-Gerlach magnets are rotated in the same axis p, = v. 

In order to connect the present discussion with section 2, we recall the definition of the spin 
trajectories as stochastic processes <f> jlig : F^ rown (i) i-> {|,|} and <^ 2t _. : F^ 2 rown (i) i-> {\,\} for times 
< t < T. We write 4>j\(2)t s — T if the j-th trajectory of the fi (r 2 ) system is in an | ensemble at 
time t. Similarly, we have 4>ji(2)t ? = i m case that the j-th trajectory of the f\ (r 2 ) system is in the [ 
ensemble at time t. 
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When the particles arrive in the space-like separated detectors at time T, they are affected by the 
forces of the magnetic fields of the Stern-Gerlach magnets. These fields depend on the axes chosen by 
the experimenters. After the particles have left the Stern-Gerlach magnets, their spin does not change 
anymore. Hence, we can describe the final result of the j'-th experiment with two axis dependent 
random variables D jip : Fj[ ovm (T) i-> {f, J.} and D j2c t : F? 2 rown (T) i-> {|, J,}, where the external forces 

F°w 2) (T) point in the same direction as the axes j2 (P). 

One gets exact anti-correlations between the spins of the j'-th particle pair during its flight if one 
can ensure for times < t < T that 

4>jit 3 + 4>j2tz- (155) 

Eq. (|155p gives rise to the entangled state in Eq. (|138() . In the viewpoint of our hidden variable theory, 
the wave-function in Eq. (|138p describes two separate ensembles, whose particles have spin states that 
are oscillating in an opposite way. 

When the entangled particles have passed the Stern-Gerlach magnets at times t > T, their spin 
values should still be anticorrelated if the axes of the devices where chosen to be the same. Accordingly, 
we have to write 

D nP # D j2p . (156) 

If a trajectory changes its ensemble from {A} ({B}) to {B} ({A}), the spin state of that trajectory 
has to change to the opposite value at the same time. Similarly, whenever a particle changes its 
spin, its {A} ({B}) ensemble must change also. Transitions like {B f} — > {A f }, {B J,} — ► {A J.}, 
{A |} — > {A J,} or {B 1*} — * {B J,} are not allowed. Otherwise, the trajectories affected would not be 
described any more by a superposition of an | and a J, state in which energy conservation holds on the 
average. 

The source of an EPR experiment repeatedly produces pairs of entangled particles. In each of these 
pairs, one particle belongs to the r\ system while the other particle is associated with vi. Therefore, 
the same index j can be used for the two trajectories of an entangled particle pair without loss of 
generality. The entangled particles start at the source with opposite spin and the same initial speed, 
due to momentum conservation. Accordingly, the trajectories that have the same index j in the f\ 
and r 2 system will be placed in the same {A} ({B}) ensemble at time t = 0, but in the opposite spin 
system. Now assume that, because of some reason, the velocities of these paired particles coincide at 
times < t < T: 



x 



• A ( B ) n\ - ^ a (b) 



iT(l)i(2)W - x ji(T)2(i)( i )- ( 157 ) 

Then, the ensemble average velocities of the rj and r 2 systems will also agree at times < t < T. The 
decision whether a trajectory x ■■rn\u2)(t) from the f\ (r 2 ) particle system changes its ensemble from 

{A} ({B}) to {B} ({A}) depends on the magnitude of its velocity x.^.L^Jt) and on the velocity 

distribution of all particles from f\ (F2) that are in the past cone of x.l 1[2) (t). The exchange mech- 
anism, as it was defined in section 4.1, works in the same way for all trajectories, no matter if they 
are in f\ or r 2 . In case Eq. (|157|) holds, it follows that the trajectories which have the same index j 
in the f\ and r 2 ensembles will switch their {A} ({B}) particle system in the same way during their 
flight. Each placement in another {A} ({B}) ensemble changes a particle's spin state to the opposite. 
Since the spin ensemble of the entangled particles was opposite from the beginning at t — 0, their spin 
trajectories will have opposite values 4>ju £ 7^ i>jit s for times < t < T. 

The velocities of all particles depend on the random forces F B 1 r ( ° 2 v f n (i), the external forces F e w 2 -,(i) 
and the trajectory exchange procedure. The exchange mechanism works in the same way for the 
particles in fi and r 2 . Furthermore, we must have 

F)f (t) = ££*(*) (158) 
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for all times < t < T if we are to observe exact correlations. Hence, any remaining differences in the 
velocity of two entangled particles can be caused by the random forces F 3 ^?^" only. Provided that we 
impose the condition 

pBrown ^ = pBrown ^ (159) 

for all j and all times t, Eq. (|15T[) and thus Eq. (|155[) hold at least for times < t < T. If the 
experimenters had chosen the same axes fi = v at the Stern-Gerlach magnets, the external forces 
induced by these devices are equal. The particles arrive in the magnets at time T. Thereby, Eq. (|155|) 
also holds for time T and all j. After the particles have passed the Stern-Gerlach magnets, their spin 
values do not change anymore and we can replace Eq. (|155p by the random variables in Eq. ()15G|) . 
We conclude that exact correlations between the spins measurement results of particles with the same 
index j in r\ and T2 can be observed. 

On the other hand if the experimenters have chosen different axes at the Stern-Gerlach magnets, 
the external forces for the j-th particles are different at time T, i. e. 

F«f (T) jk F°?(T). (160) 

In that case, the two particles which have the same index j in the f\ and r^ ensembles are subject to 
different equations of motion. It then follows from the discussion above that spin measurement results 
will not be exactly correlated for each j-th particle pair." 

In order to fulfil the requirement of Eq. (|159p . a theory must not necessarily incorporate a signalling 
mechanism between the unclupled particle systems t\ and ?2- This can be easily seen by an explicit 
construction of a local model in which Eq. (|159p holds. The only constraint on the random forces was 
that they should have a Gaussian distribution. Now consider the following mechanism: A real random 
number generator sets a value Xj € R for each j-th particle pair at preparation stage. This value Xj 
is similar to the hidden parameter in Bell's original work [5], but we define it to be independent for 
each entangled particle pair generated by the source. In our model, Xj serves as the starting value 
for two pseudo random number generators of the same type. Pseudo random number generators are 
deterministic procedures which always produce the same series provided that they are initialised with 
the same starting value Xj. After they were initialised at t = 0, we let the two generators produce 
the Gaussian distributed and time dependent series F? 1 rown (i, Xj) and F^ 2 rown (£, Xj). These series are 
the random forces of the j-th trajectories for the r\ and r-i particles. Since they only depend on their 
common starting value Xj £ R, Eq. (|159|) holds at each time. There is no signalling mechanism between 
the two disconnected forces because they simply depend on their common past. 

At first sight, one may think that our procedure would set up a deterministic passively local 
theory. Indeed, for F?-f own (£, Xj) and F^ 2 rown (i, Xj) there is a parameter Xj at preparation stage which 
determines all later outcomes of these forces. However, one does not have access to the outcomes of 
the Gaussian distributed random forces in the probability space, but only to the two spin values in 
the state space of 4>jip, ® fyjiv- 

The proof in section 3 was with regards to elements of the sigma algebras Tjs^ ^~jip and Tj2V- 
These sigma algebras contain those events on the probability space that determine the accessible spin 
results on the state space of the EPR experiment. The set of outcomes in the probability space for the 
Gaussian distributed force fields F^?^. n (t, Xj ) is considerably larger than the set of outcomes in the 
state space because the latter consists of only two spin values. Hence, every element in Tjs, Fjipi and 
Tj2V must contain many outcomes of the force fields which all lead to the same spin result in the state 
space. This implies that the condition of deterministic passive locality corresponds to the following 
situation: 

For every entangled particle that flies to a detector, an event on the probability space of the 
random force field has to happen at the first time step. All outcomes in this event must lead to 
the same particular spin state of the particle at that time. By deterministic passive locality, the spin 
result of a given particle at preparation stage has to determine the particle's spin value at the detector. 
Accordingly, the event for the Brownian force field at the first time step has to determine another event. 
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The latter consists of outcomes of the random force field at a time when the particle has arrived in 
the Stern-Gerlach magnet. Most notably, all outcomes in this predetermined event must be such that 
they lead to the same particular spin state of the particle in the detector. We will illustrate below that 
deterministic passive locality can be easily violated, even though the outcomes of the time dependent 
random force may be predetermined by a hidden parameter for each entangled particle pair. 

We let the spin trajectory ^=1,1^ of a particle, which is in a pair j = 1 and will fly to detector 
1, start in the {A |} ensemble. The magnitude of the velocity of this particle at time t ~ can 

be denoted by x^L 1 | 1 (t = 0) , and the particle's spin trajectory at that time is <j)j=is,t=o g =T- In 
the following example, we will assume that the external forces are zero until the particle arrives in a 
detector. For a given random starting value \j=i G R of the pseudo random number generators, we 
denote the outcome of the Brownian force at time t = 1 by 



uBrown [ + 
b .7=1,1 [t 



hXi 



6- 



(161) 



Since 



xf—i ii(t) is a function of F?_L° w 1 n (£, Aj=i), it may be that £1 leads to a particle velocity where 



'j=i ti(^ = 1) ^ |v(ri, t = 1)|. According to the definition of the trajectory exchange procedure from 
section 4.1, we are only allowed to replace a particle from {A} to {B} if the magnitude of its velocity 
is at least higher than |v(ri,t)|. Hence, no change of the spin trajectory will take place at t — 1 and 
we can write: 



Similarly, at time t — 2, we denote the outcome of the random force by 



F^T(*2,A j=1 ) = 6 



(162) 



where £2 maybe such that 



is allowed and we can conclude that 



xf—i u(t = 2) < |v(fi, t = 2)|. In this case, no change of the spin trajectory 



3.7=1, 1,4=0* 



0j=l,l,t=lz 



Pj'=l,l,t=2 ? - 



(163) 



There could be another pair j = 2, for which the particle 1 that flies towards detector 1 starts in the 
{A f} ensemble. The forces related to this particle may be driven by another random starting value 
Xj=2 7^ Aj=i. Without loss of generality, it is possible that we get 



te(*=l,A;= a )sil»6 



(164) 



for Xj=2 at t = 1, with 5i such that 

<f>j=2,l,t=0z — <fij=2,l,t=lg- 

Yet for t = 2, we may arrive at 



i=2,Ti 



i(t = 1) < |v(fi,i = 1)|, and thereby 



(165) 



F*i < 2 T(i = 2,A j=2 ) = <5 2 »6, 



(166) 



where 82 causes a significant change in the velocity and leads to x^ 2 T1 (t = 2) > \v(fi,t = 2)|. As a 

consequence, the particle in the {A |} ensemble may be placed in the {B J.} ensemble at time t = 2. 
Without loss of generality, we will assume that such a replacement happens for the trajectory j = 2. 
Then, x^L 2 * x (i = 2) changes to x? 2 ^(t — 2) and the the spin trajectory becomes 



»j=2,M=0j 



& J=2,M 



=1? + $3 = 



j=2,l,t=2s- 



(167) 
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That this behaviour corresponds to the violation of passive locality can be easily seen if we select 
an angle 9 = 2tt at the measurement device 1. By Eq. (|130p . the spin up state of Eq. (|163[) at t = 2 
would lead with certainty to a spin up result in detector 1. Similarly, using Eq. (|13ip . a spin down 
state can be measured with certainty at the same detector for the spin state of Eq. (|167[) at t = 2. 
These are opposite measurement results for two particles j — 1 and j = 2 that both started with 
spin up at the source and made their way to the same detector 1. Hence, there can be no event at 
preparation stage that determines a later event at the measurement devices. This corresponds to the 
failure of deterministic passive locality. 

We have to remark that even if the model with the pseudo random number generators may look 
terribly artificial, it would nevertheless be difficult to test. The deterministic functions F? L r ( ^ n (t, Xj) 
depend on Xj, which is assumed to be a real random number. Therefore, the statistics of the outcomes 
measured at time t = T at the detectors would equal the statistics of real random numbers, and this 
is actually observed in entanglement experiments. 

At first sight, it seems that a natural objection against models with a common cause Xj would be 
so-called entanglement-swapping experiments [48]. In those setups, the entangled particle pairs are 
generated by two spatially separated sources. Fortunately these experiments can be described with 
the model above, but one needs to make some minor changes: The area of the particle source in section 
2.1 has to be replaced by a slice in which the two particle sources are. That slice must be disjoint from 
the future cones of the two Stern-Gerlach magnets and all proofs in section 2 and 3 work as before. 
In section 6.2, Eq. (|159|) has to hold for the Brownian forces at the two separated particle sources. 
Nevertheless, it is possible to find a spatial region where the past cones of the two sources intersect. 
At this place, a common cause Xj € R could well be established, leading to Eq. (|159|) . Hence, even 
with entangled particles that are produced by two independent sources, there is no need for hidden 
variable theories to include an instantaneous connection between spatially separated regions. 

It is well known that entanglement phenomena arc not restricted to spin states. For example, 
according to quantum mechanics, exact correlations between positions of entangled particles should be 
observed when they are emitted by the source at almost the same time. A first outline of an experiment 
to investigate this behaviour was presented in 1989 [49], and the analogue of this setup for photons was 
realised in 2005 [50]. Recently a proposal for a similar experiment with material particles has been 
made [51], which still waits for its realisation. Here, we only mention that the correlations predicted 
by quantum mechanics for these situations are naturally explained by our local model. Obviously, the 
discussion above is valid as long as the entangled particles in each pair are ejected within an interval 
smaller than the course grained timescale r = 2k h ' T of the Brownian motion process. In turn, we are 
lead to Eq. (|157[) , which describes particles that behave exactly correlated in velocity and time. After 
an integration of the particle velocities, we get exact time and position correlations. 

In order to violate Bell's inequality with a local model of the EPR experiment, it does suffice to 
assume that the Gaussian distributed random forces F? 1 rown (£) and F? 2 rown ^ have a kind of "memory" 
of their past. Nevertheless, the author has to admit the somewhat incomplete nature of this article. 
The model given above is of course terribly artificial. It serves only as an example to show that 
a local hidden variable theory which violates Bell's inequality can be build in principle. To give a 
more physical reason for Eq. (|159[) and to replace the rather arbitrary algorithm for the trajectory 
exchange procedure, one would need further studies. In fact, passive locality may be violated in 
many other local hidden variable theories, which are all sharing the property that the events at one 
measurement apparatus do not influence the events at another spatially separated device. Besides 
from Nelson, who gave two examples of such theories, a different model with similar locality features 
was published recently [52] . The authors of this article did not study Bell's theorem in the language 
of advanced mathematical probability, as it was done by Nelson, but they seem to share his opinion 
about the meaning of passive locality. Here, we wish to leave it at that with a citation by Nelson from 
a popular talk. It shows that behaviour which is similar to entanglement is in fact commonly observed 
in stochastically evolving systems, and not at all confined to quantum mechanical objects: 

"The snowflake problem is this: The inexhaustible variety of snowflakes makes it evident 
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that chance plays a major role in their development, yet they always preserve hexagonal 
symmetry - how does a portion of the snowflakc growing at random on one side know to 
grow in precisely the same fashion as its partner all the way over on the other side? This 
is mysterious but not beyond understanding" [53]. 

"Similar phenomena arc well known to occur in the study of random fields. [...] Perhaps 
the results of correlated spin experiments are ultimately no more and no less, mysterious 
than is the random growth of snowflakes with hexagonal symmetry" [9] . 

An explicit example of classical Brownian particles that exhibit entanglement- like phenomena is studied 
in [54]. 

7 Objections against Nelson's theorem 

7.1 The argument of Greenberger et al. — Bell's theorem without inequal- 
ities 

The essential argument in Nelson's (and Jarret's) reasoning is that a family of device dependent product 
random variables 4>jip <S> 4>j2H may generate events at spatially separated locations for which passive 
locality fails, i. e. 

p 0jim«"*j2p \ a jip l) a 3 2 " ^ s ) ^ p </w»<A 3 2,? Wj\p,\fjs) P 3M «.0 3 2,? {°j2s\ Fjs) ■ (168) 

Nelson and Jarret explicitly wrote that this can happen without any signalling mechanism in the 
theory that produces the events in question. The latter statement, although mathematically correct, 
has raised criticisms from physicists since it was made. For example, in his famous introductory 
textbook, Ballcntine [55] writes: 

"There have been some debates as to whether the factorisation of the probability is justified by 
the locality condition alone, or whether it requires some additional stronger assumptions. If so, the 
contradiction may be blamed on those additional assumptions. The issue is subtle, but fortunately it 
is now irrelevant since the new proof in [56] does not make use of probability" . 

The assertion that the work [56] by Greenberger et al. does not make use of probability is true. 
However, their arguments do only account for deterministic hidden variable theories of a sort where 
the events at the particle source are equivalent to the ones at the detectors. The authors are aware of 
this, noting at page 1138 in their article: 

"The salvaging strategy could consist in [...] a stochastic local theory in contrast to a deterministic 
local theory which has been considered so far in this paper" . 

They then go on, writing that these "stochastic local theories" would be impossible. They claim that 
this would have been thoroughly investigated by Bell [57], and Clauscr and Home [58]. Furthermore, 
they mention so-called "equivalence theorems" , which were proven by Stapp [59] and Fine [60] . Those 
theorems are said to show that each stochastic local hidden variable theory could be duplicated by 
a deterministic passively local one. However, all articles cited by [56] only consider theories where 
passive locality is assumed to hold. (Bell [57] uses passive locality explicitly in his Eq. (5), Clauser and 
Home [58] make this assumption in their Eqs. (2) and (2'), Stapp [59] does the same in his Eq. (7a) and 
also Fine [60] on pp. 1306 and 1308.) In fact, an equivalence to a deterministic passively local theory 
can only be proven in this case. We conclude that the proof in [56] does not call actively local hidden 
variable theories into question in which passive locality fails. That the arguments of Greenberger et 
al. do not apply for stochastic hidden variable theories has also been shown by Khrennikov [61] in a 
detailed analysis. 
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7.2 The argument of Diirr — Bell's theorem with a single probability mea- 
sure P 

Furthermore, there are researchers who try to derive Bell's theorem in the form of Eq. (|37p . but without 
the assumption of passive locality. They (e. g. [62] on p. 182) do not write the inequality with a whole 
family of axis dependent probability measures P^, jliX ®$ j2P for device dependent events (Jjifi €E Fjifi and 
<Jj2v € Fjiv- Instead, they are using one single detector independent measure P for events that are 
thought to be produced by the source and, nevertheless, as dependent on the axes of the detectors. 
This is problematic because if we compute the quantum mechanical correlation coefficient 

EC/r^s^K/w)!®^!*) (169) 

with the entangled state in Eq. (|138[) . we are definitely led to Eq. (fl6|) and to a family of axis dependent 
measures P^ 1?( g><^ 2i7 . 

By using one probability measure for axis dependent events generated by the particle source, we 
would only be able to investigate theories where an intrinsic particle property is defined at preparation 
stage that determines the later results at the two detectors for each conceivable pair of axes. Yet it 
may be possible that the source generates particles whose properties are modified later by the Stern- 
Gcrlach magnets. In this case, a hidden variable model of the EPR experiment has to use an enlarged 
probability space. The events at preparation stage take place in a smaller space and their probabilities 
do not depend on the detector's axes, whereas the events at the detectors happen in the enlarged 
probability space. A natural assumption is that the events at the two separated measurement stations 
are generated by a family of device dependent product random variables 4>j\p, <£> 4>j2P- Accordingly, 
we arrive at a family of axis dependent probability measures P^ 1 3®^ J2P for the joint probabilities of 
the events at the two Stern-Gerlach magnets. Since the probability measures from the formalism of 
quantum mechanics are of the same structure, we believe that Nelson's mathematical analysis of the 
EPR experiment is the correct one. 

At this point, we want to remark that the impossibility to describe the EPR experiment with one 
single probability space, and the need for the inclusion of the detector settings in the corresponding 
random variables was also mentioned recently by Khrennikov [63] . 

Nevertheless, it is possible to confine oneself to a single and axis independent measure P, as long 
as active and passive locality are explicitly stated to hold. With that proposition and the proof in 
section 3, one can trace the events <Tj\pt € fjijl an( i a j2p, € Fj2p, f° r an ax is P at the detectors back 
to equivalent events in Tjs at the common preparation stage. We can repeat this procedure with two 
other pairs of equivalent events for different axes v and V . Eventually, after some further calculations, 
we will arrive at an inequality which is analogue to Eq. (|37p but defined only for events in Tjg- Due 
to active locality, the events at preparation stage can not be modified later by the settings of the 
detectors. Hence, with the assumptions of active and passive locality, we can use an axis independent 
probability measure P to analyse the EPR experiment. A proof of Bell's theorem that goes exactly 
this way can be found in [12]. 

7.3 The argument of Redhead - stability against disturbances 

A different argument against stochastic hidden variable theories was raised by Redhead [26] (who 
presents in his book at p. 98 a proof of Bell's inequality which is very similar to that of Nelson, 
although with some differences). Redhead argues that stochastic hidden parameter theories would be 
impossible and advances this opinion with the following observation: He asserts that under natural 
conditions, one can never exclude small disturbances, which might happen erroneously to a single- 
particle (maybe through a random disturbance by the gases of the air in the laboratory). Redhead 
then defines what he calls "the most general perturbation that is not a phase shift" as 

u(n^) = e (sn ^/ 2 , (170) 
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Figure 5: Velocity distribution of the particles in the {A} ensemble. Only particles whose velocity is 
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8ji(t) around |v(ri, i)\ will be interchanged only because of the 



where n is an arbitrary unit vector and d is an angle about < $ < Att. The perturbation should be 
confined to one part of the overall Hilbert space, e. g. H.2 only. Accordingly, the entangled state in 
Eq. (|138p will transform to 



y'(r u r 2 ,t) ^—(Mri,t)\-) ®u{n 7 d)Mr2,t)\+)) 



-1 



+ -= (Mn, t) |+) <g> u(n, tf)V>i(r 2 , t) |-» . 



(171) 



Redhead concludes that the exact correlations observed in EPR experiments will be lost with 
^'(ri, r*2,£) for certain axes of the detectors. Nevertheless, the phenomenon of entanglement can be 
experimentally verified through distances over 100 kilometres in the air [64], where small perturbations 
acting on single particles cannot be avoided. 

Yet, in spirit of Ballentine's ensemble interpretation [31], a quantum mechanical state may not 
describe individual particles, but rather an infinitely large statistical ensemble of trajectories. In this 
case, a perturbation like Eq. (|170p would correspond to a disturbance of an entire statistical ensemble 
of infinitely many particles for a sufficiently long time. Moreover, if defined as above, u(n, $) acts on 
every particle of that ensemble in exactly the same way. This has nothing to do with a small erroneous 
disturbance that may happen accidentally to individual particles in the laboratory. 

In this context, we want to show that our theory is, in some sense, especially robust against small 
disturbances, at least if they act on individual particles. For a long distance entanglement experiment, 
it is a reasonable requirement to expect large uncorrelated external forces F2 Xt (r*2,i) and Ff ct (ri,t), 
which act on the particles during their flight, to be absent. It is well known that without external 
forces, the Langevin equations ([55)) lead to a Maxwell distribution of the particle's velocity. Hence, 
most fi(2) particles in the {A} or {B} ensembles have a velocity around one centre w v^ (^1(2), t) or 

~|v?(r 1(2) ,i)|- 

Without loss of generality, we will assume that a small disturbance Sji(t) induces a slightly different 
velocity of the j-th f\ particle in the {A |} ensemble at a point (fi,t) in space-time: 



&(*) 



= s£ n (t) 



5 3l (t). 



(172) 



We define 6ji(i) to be small if the half-width of the particle's velocity distribution is by far larger 



than 



M*) 



Moreover, the disturbance should be small enough that the perturbation can't change 



the statistics of the majority of particles through the trajectory exchange mechanism. The latter 



:-!(» 



requirement is equivalent to the condition 
||v(f 1 ,<)|-|vA(r 1 ,t)||» 5ji(t) . 



(173) 



Additionally, we expect the components of Sj±(t) to be fluctuating with zero mean. 

A particle can only be replaced from the {A f} to the {B J.} ensemble if the magnitude of its 

velocity xA (£) is higher than |v(fi,£)|. The disturbance does not contribute to |v(r*i,£)| because the 
average of each component in Sj\(t) is assumed to be zero. Hence, 



\y(n,t)\- 4 n (t) < 5 3l {t) 



(174) 



has to hold if the perturbation Sj\(t) is the primary reason for x^-^i) to be placed in an ensemble. 
The particles fulfilling Eq. (|174|l are in a small velocity interval of width 5ji(t) around |v(ri,t)| (see 

Fig. ED. 



Since 



M*) 



is assumed to be much smaller than the half-width of the particle's velocity distri- 
bution, the number of trajectories for which Eq. (|174[) holds will be low. Due to equation (|173[) . the 
interval in which the disturbance is able to change the ensemble of a particle is far away from the 
centre |v^(fi,£)| where most trajectories are found. This further reduces the number of particles that 

can be affected by Sji(t). 

All in all, we can conclude that the number of particles whose ensemble is changed solely because 
of the disturbance Sji(t) is very low. For this reason, the disturbance can only cause a vanishingly 
small reduction of the efficiency of the entanglement in an EPR experiment, but not its breakdown. 
Hence, the theory we have constructed is especially robust against small disturbances unavoidable in 
a laboratory, in contrast to the claims by [26] . 



8 Conclusions and outlook 

8.1 What has been achieved 

The most popular argument against hidden variable theories with particle trajectories is, as commonly 
taken from Bell's words [that those theories] "require a mechanism, whereby the setting of one mea- 
suring device can influence the reading of another device" and "moreover, the signal involved must 
propagate instantaneously" . As shown by Nelson, whose analysis we reviewed here, this assertion is 
physically unfounded. Specially designed stochastic processes can violate Bell's inequality at two sep- 
arated locations without any contact at all. We have mathematically analysed, what such a behaviour 
really does imply for a stochastic theory, and we have constructed a corresponding one. 

8.2 Possibility of relativistic extension 

Unfortunately, our theory has, at least to this day, no relativistic extension. Stochastic mechanics in 
the original form in which Nelson wrote it down had severe conceptual difficulties in the relativistic 
domain. These problems arose due to the Markov property of the stochastic process one had to use. 
In fact, it was shown that the definition of a Markov process on Minkowski space would lead either 
to determinism or to the violation of causality [65-67]. Nevertheless, with some different assumptions 
for the stochastic process involved, the Klein-Gordon equation in four dimensions [68-73] as well as 
the Dirac equation in 1+1 dimensions [74-76] could be derived in various investigations. Furthermore, 
there also exist attempts to derive the Dirac equation in 3 + 1 dimensions from various versions of 
stochastic mechanics e. g. [77,78]. Yet, the derivation of the Dirac equation from an actively local 
stochastic process still seems to be a completely open problem. Therefore, we will try to analyse if 
some of the methods used by Fritsche and Haugk could perhaps help to solve it. 
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At first, we want to emphasise that unlike Nelson's original formulation, the theory of Fritsche 
and Haugk describes a non-markovian process and it thereby avoids the no-go theorems of [65-67]. 
However, in order to rewrite the Fritsche-Haugk model in a manifestly covariant form, one would need 
a fully formulated theory of special relativistic Brownian motion. Regrettably, although such a theory 
certainly must exist (since non-relativistic Brownian motion does), it seems that till this day, no one 
gave a suitable model based on a relativistic generalisation of the Langevin equations. (Although there 
are highly interesting attempts in this direction [79-82].) 

Fortunately, there are other ways to describe a relativistic system that evolves stochastically. It 
is shown in [15] on p. 386 that after an approximation, the Eqs. in (|56[) can be transformed into the 
Navier-Stokes equations. We have in Eq. (|56| for the {B} system: 

I- (v B - u B ) + (v B + u B ) • V (v B - it B ) - uA (v B - u B ) = — F oxt (175) 

Ol 772q 

If we set Vp B (r, t) = 0, we observe that 



u 



vV\n(p B (r,t)/p ) = -- — — Vp B (f,*) = 0. (176) 



With p(r,t) denoting the pressure and p m = mop B (r, t) as the mass density, the following relation 
holds: 

p m F 0Xt = -Vp(r,t). (177) 

Inserting the result of Eq. (|176[) and the relation (|177p into Eq. (|175[) . we arrive at the Navier-Stokes 
equation of a fluid 

Pm^v B + Pm (v B • Vv B ) - J/p m Av B + Vp(f, t) = 0. (178) 

For Eq. (|178p . various manifestly covariant formulations have been proposed, see [83-86], and for 
a general review [87]. In order to derive the relativistic version of Eq. (|178|) from first principles, one 
typically starts with relativistic Boltzmann equations [88-95] . Due to the close similarity between Eq. 
IJ178JI and Eq. (|56|) . it seems at least plausible that it might be possible to formulate the theory of 
Fritsche and Haugk in a relativistic setting. 
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